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SUMMARY 

This  study  performs  maximum  likelihood  event  magnitudes  and  station  effects  for 
111  nuclear  explosions  from  more  than  8  different  test  sites  using  a  general  linear 
model.  The  estimates  are  performed  for  the  three  P-wave  phases  "a",  "b",  and  "max" 
by  considering  the  signals,  non-detection,  and  clipping  information.  The  bias  between 
a  model  that  does  not  include  censoring  effects  from  non-detection  and  clipping  is 
demonstrated.  The  maximum-likelihood  magnitude  estimates  are  performed  using  data 
for  each  phase  individually  and  the  three  phases  combined.  The  difference  in  mj,’s 
between  the  phases  is  studied  for  some  of  the  major  test  sites.  A  presumed  multiple 
event  is  found  to  be  characterized  with  large  "max"  phase  amplitude  compared  to  the 
"a"  phase,  whereas  in  the  case  of  a  presumed  cratering  event  showing  surficial  col¬ 
lapse,  the  amplitudes  for  the  "b"  and  "max"  phases  are  smaller  instead.  In  addition, 
we  also  examined  the  biases  at  several  test  sites  that  may  exist  for  general  linear 
models  based  on  a  specific  test  site  relative  to  a  model  that  contains  sources  from 
many  test  sites.  The  effects  of  assuming  correlation  among  the  events  on  the  magni¬ 


tude  estimates  using  the  present  event-station  distribution  are  investigated  using  an 


inter-event  correlation  model. 
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mi,(a)  event  magnitudes  from  LS-GLM  (no-censoring)  estimates  24 
plotted  against  MLE-GLM  (censoring)  estimates  for  106  explo¬ 
sions.  Note  that  there  is  only  a  narrow  range  of  magnitudes  (5.7 
to  6.5)  for  which  the  two  sets  of  magnitudes  are  on  average  the 
same.  The  LS-GLM  event  magnitudes  are  too  large  at  the  lower 

mb’s. 

mb(b)  event  magnitudes  from  LS-GLM  (no-censoring)  estimates  25 
plotted  against  MLE-GLM  (censoring)  estimates  for  111  explo¬ 
sions.  Note  that  there  is  only  a  narrow  range  of  magnitudes  (5.5 
to  6.5)  for  which  the  two  sets  of  magnitudes  are  on  average  the 
same.  The  LS-GLM  event  magnitudes  are  too  large  at  the  lower 
mb’s  and  too  small  at  the  higher  mb’s. 

mb(max)  event  magnitudes  from  LS-GLM  (no-censoring)  esti-  26 
mates  plotted  against  MLE-GLM  (censoring)  estimates  for  111 
explosions.  Note  that  there  is  only  a  narrow  range  of  magnitudes 
(5.5  to  6.5)  for  which  the  two  sets  of  magnitudes  are  on  average 
the  same.  The  LS-GLM  event  magnitudes  are  too  large  at  the 
lower  mb’s  and  too  small  at  the  higher  mb’s. 

LS-GLM  (no-censoring)  mb’s  estimates  for  the  "a"  phase  plotted  28 
against  the  "b"  phase.  The  mb(b)'s  are  on  the  average  0.16  mag¬ 
nitude  unit  greater  than  the  mb(a)’s  for  the  larger  events. 

LS-GLM  (no-censoring)  mb’s  estimates  for  the  "a"  phase  plotted  29 
against  the  "max"  phase.  The  mb(max)’s  arc  on  the  average  0.4 
magnitude  unit  greater  than  the  mb(a)’s  for  the  larger  events. 

LS-GLM  (no-censoring)  mb’s  estimates  for  the  "b"  phase  plotted  30 
against  the  "max"  phase.  The  mb(max)’s  are  on  the  average  0.2 
magnitude  unit  greater  than  the  mb(a)’s  for  the  larger  events. 

MLE-GLM  mb’s  estimates  for  the  "a"  phase  plotted  against  those  31 
for  the  "b"  phase.  These  mb’s  are  derived  using  aU  three  phases. 

The  mb(b)’s  arc  on  the  average  0.28  magnitude  unit  greater  than 
the  mb(a)’s. 
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3b  MLE-GLM  mb’s  estimates  for  the  "a"  phase  plotted  against  those  32 

for  the  "max"  phase.  These  mb’s  are  derived  using  all  three 
phases.  The  mb(max)’s  are  on  the  average  0.48  magnitude  units 
greater  than  the  mb(a)’s. 

3c  MLE-GLM  mb’s  estimates  for  the  "b"  phase  plotted  against  those  33 

for  the  "max"  phase.  These  mb’s  are  derived  using  all  three 
phases.  The  mb(max)’s  are  on  the  average  0.19  magnitude  units 
greater  than  the  mb(a)’s. 

4a  LS-GLM  mb(a)  estimates  derived  from  using  just  the  "a"  phase  35 

plotted  versus  those  derived  from  using  all  tliree  phases.  The 
mb(a)’s  derived  from  using  just  the  "a"  phase  are  consistently 
higher. 

4b  LS-GLM  mb(b)  estimates  derived  from  using  just  the  "b"  phase  36 

plotted. versus  those  derived  from  using  all  three  phases.  These 
two  sets  of  magnitudes  agree  very  well. 

4c  LS-GLM  mb(max)  estimates  derived  from  using  just  the  "max"  37 

phase  plotted  versus  those  derived  from  using  all  three  phases. 

These  two  sets  of  magnitudes  agree  very  well. 

5a  MLE-GLM  mb(a)  estimates  derived  from  using  just  the  "a"  phase  38 

plotted  versus  those  derived  from  using  all  three  phases.  There  is 
a  slight  indication  that  the  mb(a)’s  for  just  the  "a"  phase  are 
biased  low  towards  low  mb’s. 

5b  MLE-GLM  mb(a)  estimates  derived  from  using  just  the  "b"  39 

phase  plotted  versus  those  derived  from  using  all  three  phases. 

These  two  sets  of  magnitudes  agree  very  well. 

5c  MLE-GLM  mb(max)  estimates  derived  from  using  just  the  "max"  40 

phase  plotted  versus  those  derived  from  using  all  three  phases. 

These  two  sets  of  magnitudes  agree  very  well. 

6a  MLE-GLM  station  effects  derived  using  the  "a"  phase  plotted  46 

versus  those  for  the  "b"  phase.  The  bias  between  the  two  esti¬ 
mates  is  up  to  0.2  magnitude  units. 

6b  MLE-GLM  station  effects  derived  using  the  "a"  phase  plotted  47 

versus  those  for  the  "max"  phase.  The  bias  between  the  two  esti¬ 
mations  is  much  higher,  being  up  to  0.4  magnitude  units. 
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6c  MLE-GLM  station  effects  derived  using  the  "b"  phase  plotted  48 

versus  those  for  the  "max"  phase.  The  bias  between  the  two  esti¬ 
mates  is  up  to  0.2  magnitude  units. 

7  LS-GLM  station  effects  versus  MLE-GLM  station  effects  for  127  49 

WWSSN  stations.  Note  that  for  most  stations  the  difference  is 
skewed  and  therefore  potential  for  bias  exists  when  the  censoring 
information  is  ignored. 

8  MLE-GLM  mf,  error  estimates  for  using  bootstrap  versus  EM.  52 

The  bootstrap  error  estimates  are  biased  high. 

9a  LS-GLM  mb  error  estimates  versus  MLE-GLM  nib  error  esti-  53 

mates  obtained  from  EM.  The  LS-GLM  mb  errors  computed  with 
no-censoring  information  are  higher  than  those  with  censoring  in¬ 
formation. 

9b  LS-GLM  mb  error  estimates  versus  MLE-GLM  nib 

mates  obtained  using  bootstrap.  The  LS-GLM  mb  error  estimates 
computed  with  no-censoring  information  are  higher  than  those 
with  censoring  information. 

10a  LS-GLM  nib  error  estimates  versus  MLE-GLM  nib  ®sti-  56 

mates  obtained  from  EM  using  just  the  "a”  phase  data.  The  LS- 
GLM  mb  error  estimates  computed  with  no-censoring  informa¬ 
tion  are  higher  than  those  with  censoring  information. 

10b  LS-GLM  mb  error  estimates  versus  MLE-GLM  nib 

mates  obtained  from  EM  using  just  the  "b"  phase  data.  The  LS- 
GLM  mb  error  estimates  computed  with  no-censoring  informa¬ 
tion  are  higher  than  those  with  censoring  information  but  show 
smaller  bias  than  those  for  phase  "a". 

10c  LS-GLM  rrib  error  estimates  versus  MLE-GLM  nib  error  esti-  58 

mates  obtained  from  EM  using  just  the  "max"  phase  data.  The 
LS-GLM  mb  error  estimates  computed  with  no-censoring  infor¬ 
mation  are  higher  than  those  with  censoring  information  but 
show  smaller  bias  than  those  for  phase  "a". 

11a  LS-GLM  site  term  error  estimates  ve'sus  MLE-GLM  site  term  59 

error  estimates  using  EM.  The  error  estimates  for  the  no- 
censoring  case  are  larger  than  those  for  the  censoring  case  in 
agreement  with  the  mb  error  observations. 
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MLE-GLM  site  term  error  estimates  obtained  using  bootstrap  are  60 
plotted  against  those  obtained  using  EM.  The  error  estimates  us¬ 
ing  the  two  different  schemes  are  in  quite  close  agreement  to 
each  other. 

mb(max)’s  plotted  against  Smt’s  for  Shagan.  The  m^’s  are  max-  63 
imum  likelihood  estimates  with  the  censoring  information.  The 
solid  circles  are  mb(max)-m5(a).  The  open  circles  are  mb(max)- 
mb(b). 

mb(max)’s  plotted  against  Smb’s  for  Degelen.  The  mb’s  are  max-  64 
imum  likelihood  estimates  with  the  censoring  information.  The 
solid  circles  are  mb(max)-mb(a).  The  open  circles  are  mb(max)- 
mb(b).  The  hatched  circles  arc  mb(max)-mb(a)  for  which  there 
are  less  than  ten  "a"  phase  signal  readings. 

mb(max)’s  plotted  against  Smb’s  for  NTS.  The  mb’s  are  max-  66 
imum  likelihood  estimates  with  the  censoring  information.  The 
solid  circles  are  mb(max)-mb(a).  The  open  circles  arc  mb(max)- 
mb(b).  The  hatched  circles  are  mb(max)-mb(a)  for  which  there 
are  less  than  ten  "a"  phase  signal  readings. 

mb(max)’s  plotted  against  Smb’s  for  Novaya  Zemlya.  The  mb’s  67 
are  maximum  likelihood  estimates  with  the  censoring  informa¬ 
tion.  The  solid  and  open  circles  are  mb(max)-mb(a)  and 
mb(max)-mb(b),  respectively,  for  the  northern  site.  The  solid  and 
open  squares  are  mb(max)-mb(a)  and  mb(max)-mb(b),  respective¬ 
ly,  for  the  southern  site.  The  hatched  circles  and  squares  arc 
mb(max)-mb(a)  for  which  there  are  less  than  ten  "a"  phase  signal 
readings. 

Event  magnitudes  for  phase  "a",  mb(a),  from  MLE-GLM  solu-  74 
lions  assuming  event-event  correlation  versus  no  event-event 
correlation. 

Event  magnitudes  for  phase  "b",  mb(b),  from  MLE-GLM  solu-  75 
tions  assuming  event-event  correlation  versus  no  event-event 
correlation. 

Event  magnitudes  for  phase  "max",  mb(max),  from  MLE-GLM  76 
solutions  assuming  event-event  correlation  versus  no  event-event 
correlation. 
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nib’s  estimates  assuming  event-event  correlation  for  the  "a"  phase  78 
plotted  against  the  "b"  phase.  The  bias  is  similar  to  that  observed 
for  the  MLE-GLM  case. 

mb’s  estimates  assuming  event-event  correlation  for  the  "a"  phase  79 
plotted  against  the  "max"  phase.  The  bias  is  similar  to  that  ob¬ 
served  for  the  MLE-GLM  case. 

mb’s  estimates  assuming  event-event  correlation  for  the  "b"  80 

phase  plotted  against  the  "max"  phase.  The  bias  is  similar  to  that 
observed  for  the  MLE-GLM  case. 

Site  effects  for  MLE-GLM  solutions  assuming  event-event  corre-  8 1 

lation  versus  LS-GLM  with  no-censoring  information  based  on 
mb’s  for  alt  phases. 

Site  effects  for  MLE-GLM  solutions  assuming  event-event  corre-  82 
lation  versus  no  event-event  correlation  based  on  mb’s  for  all 
phases. 
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1.0  INTRODUCTION 

This  study  presents  results  of  the  analysis  of  the  P-wave  amplitudes  for  111 
explosions  from  over  8  major  test  sites  world-wide  as  well  as  numerous  US  and 
foreign  peaceful  nuclear  explosions  located  away  from  the  established  test  sites.  The 
data  base  consists  of  seismic  measurements  made  from  the  short-period  WWSSN  film 
chips  for  these  events.  The  maximum-likelihood  estimation  general  linear  model  pro¬ 
gram,  MLGLM,  was  developed  to  estimate  the  magnitudes  for  these  events  in  a  model 
framework.  In  magnitude-yield  estimation  studies,  it  is  important  to  provide  event- 
magnitude  calibration  within  a  large  data  set  to  give  relative  magnitudes  between 
events.  We  report  on  the  estimation  of  maximum-likelihood  magnitudes  and  station 
correction  terms,  MLE-GLM88.  This  model  is  used  to  document  some  of  the  statistical 
biases  that  occur  in  typical  processing  of  data  for  magnitude:yield  estimation.  In  addi¬ 
tion,  we  have  also  incorporated  event-event  correlation  into  the  magnitude  estimations 
for  the  large  database,  which  involves  inverting  an  immense  matrix  with  the  aid  of  a 
supercomputer.  A  maximum-likelihood  estimation  of  the  error  terms  is  also  used  which 
proves  to  be  more  efficient  and  practical  than  the  bootstrap  method. 

Douglas  (1966)  has  applied  the  general  linear  model  (GLM)  to  magnitude  estima¬ 
tion  as  the  least  squares  multiple  factor  program  (LSMF).  Heasler  et  al.  (1988)  have 
discussed  the  formulation  of  a  modeling  and  estimation  theory  (GLM)  and  tested  it  on 
NTS  data.  The  GLM  program  is  widely  used  to  calibrate  a  network  of  stations  and 
gives  estimates  of  the  precision  of  relative  event  magnitudes  for  a  given  test  site.  The 
model  assumes  that  the  observed  magnitude  at  each  station  for  each  event  is  the  sum 
of  independent  event  and  station  terms: 
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mjj  —  Ej  +  Sj  +  Ejj 

where  Ej  is  the  i’th  event  magnitude,  Sj  is  the  j’th  station  term,  and  Cjj  is  an  indepen¬ 
dent  normal  random  error  with  zero  mean  and  standard  deviation,  a.  This  linear 
model  has  been  validated  and  solved  using  standard  linear  methods.  However,  tradi¬ 
tional  solutions  to  this  system  of  linear  equations  cannot  accommodate  censoring  infor¬ 
mation  such  as  clipping  or  non-detection  levels  for  the  individual  station  magnitudes. 
Such  information  has  been  found  to  be  relevant  to  the  magnitude  estimation  problem 
either  when  events  are  small  and  low  amplitude  stations  are  excluded  by  non -detection, 
or  when  events  are  large  and  large  magnitudes  are  excluded  by  clipping.  In  either 
case,  we  can  use  the  censoring  information  such  as  mjj  <  m'^jj,  or  mjj  >  m'^jj,  in  a  max¬ 
imum  likelihood  formulation  of  the  GLM,  where  m‘*jj  is  the  magnitude  at  the  lowest 
detectable  level  and  m®  j  is  the  magnitude  at  the  critically  clipped  level  to  provide  con¬ 
straints  to  the  model. 

Ringdal  (1977)  introduced  the  maximum  likelihood  estimator  (MLE)  to  correct 
for  the  statistical  bias  introduced  by  data  censoring  from  non-detection.  Von  Seggem 
and  Rivers  (1978)  extended  the  procedure  to  include  data  censoring  by  clipping. 
Blandford  and  Shumway  (1982)  further  extended  the  maximum  likelihood  estimator  to 
the  GLM  in  the  presence  of  data  censoring  using  the  expectation  maximization  algo¬ 
rithm  (MLE-GLM).  They  simultaneously  estimated  event  magnitudes  and  station 
corrections  in  a  maximum  likelihood  sense.  This  MLE-GLM  is  to  the  LSMF  program 
what  the  Ringdal  MLE  procedure  is  to  the  single  event  network  average.  The  follow¬ 
ing  table  summarizes  the  inter-relationships  of  these  estimation  procedures.  In  the  case 
of  no  censoring,  the  MLE-GLM  reduces  to  the  least  squares  linear  estimation  problem 
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(LS-GLM)  and  in  the  case  of  no  station  effects  the  least  squares  linear  problem 
reduces  to  a  network  average. 


Table  1.  Inter-Relationships  of  Magnitude  Estimators 

No  station  effects 

1  Dimension 

station  effects 
N+M-1  Dimensions 

No  censoring 

Network  Mean 

LS-GLM 

Censoring 

MLE 

MLE-GLM 

where  M  is  the  number  of  stations  and  N  is  the  number  of  events. 

The  expectation  maximization  algorithm  was  applied  to  the  linear  estimation 
problem  with  two  additional  improvements  by  Blandford  et  al.  (1984),  They  proposed 
to  use  a  bootstrapping  scheme  to  estimate  the  uncertainty  for  the  maximum  likelihood 
general  linear  model  and  to  introduce  data  correlation  estimates  to  the  GLM  in  order  to 
reduce  the  biasing  effects  of  uneven  data  sampling. 

The  general  linear  model  (GLM)  was  improved  by  Blandford  and  Shumway 
(1982)  to  relate  amplitude  readings,  event  magnitudes,  and  station  effects  to  include 
the  data  censoring  effects  of  clipping  and  non-detection.  The  technique  was  applied  to 
WWSSN  short-period  recordings  of  NTS  and  Algerian  shots  in  granite  (Blandford  and 
Shumway,  1982),  and  after  conections  for  pP  it  was  concluded  that  there  was  only  a 
small  offset  in  magnitude  between  the  NTS  and  Algerian  shots.  A  number  of  other 
event  magnitudes  were  also  determined,  and  the  station  corrections  resulting  were 
compared  to  those  of  North  (1977)  with  fairly  good  agreement. 
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It  was  noted  in  Blandford  and  Shumway  (1982)  that  a  weakness  of  the  mathemat¬ 
ical  procedure  was  the  assumption  that  readings  from  explosions  in  close  proximity  to 
each  other  were  independent.  Thus  if  many  events  were  entered  into  the  linear  model 
from  a  single  test  site,  and  if  those  amplitudes  were  in  fact  highly  correlated,  then  the 
station  corrections  would  tend  to  reflect  that  particular  path  and  not  truly  represent  the 
average  structure  under  the  stations.  Again,  this  could  lead  to  inaccurate  relative  mag¬ 
nitudes  between  events  at  widely  different  test  sites.  In  Blandford  et  al.  (1984)  this 
was  remedied  by  estimating  the  event-event  correlation  and  incorporating  the  resulting 
correlation  matrix  (regarded  as  a  priori )  into  the  theory. 

A  method  of  computing  confidence  intervals  for  this  complex  problem  was  also 
developed  and  applied  by  Blandford  et  al.  (1984)  using  the  bootstrap  methods  of 
Efron  (1979,  1981).  This  becomes  an  attractive  method  for  uncertainty  estimation  as 
computing  power  becomes  less  expensive  and  more  readily  available.  It  is  an  extension 
of  the  "Jackknife"  method  where  each  N-1  subset  of  the  N  data  is  examined  in  turn  to 
test  for  stability  of  the  results  to  subtle  changes  in  the  data  set.  The  bootstrap  pro¬ 
cedure  resamples  the  entire  data  set  M  times  treating  the  data  as  the  next  best  thing  to 
an  ensemble  of  realizations.  Each  resampling  is  used  to  estimate  the  model  parameters 
and  the  ensemble  of  model  parameters  is  used  to  estimate  the  stability  of  the  final 
model  estimate.  This  procedure  is  used  to  evaluate  the  results  of  the  maximum  likeli¬ 
hood  estimators  for  the  GLM  rather  than  evaluate  the  multi-dimensional  partial  deriva¬ 
tive  matrix  of  the  GLM’s  likelihood  function.  McLaughlin  (1988)  has  shown  how  to 
apply  the  bootstrap  procedure  to  single  event  maximum  likelihood  magnitude  estima¬ 
tion.  Of  course,  the  bootstrap  procedure  has  many  more  potential  geophysical 
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applications  than  just  the  magnitude  estimation  problem. 
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2.0  DATA  ANALYSIS 

The  expectation  maximization  (EM)  algorithm  (Dempster  et  al. ,  1977)  may  be 
applied  to  solve  the  multi-parameter  version  of  the  maximum  likelihood  estimation 
problem  originally  considered  by  Ringdal  (1977)  when  the  errors  in  observing  the 
amplitudes  arc  uncorrelated  from  event  to  event  and  from  station  to  station.  This 
method,  as  has  been  applied  by  Blijidford  and  Shumway  (1982)  is  desirable  since  the 
log-likelihood  is  a  non-linear  function  of  the  residuals,  non-detection  levels,  and  clip¬ 
ping  levels.  The  number  of  unknowns  (N+M-1)  for  the  multi -parameter  estimation 
problem  is  equal  to  one  less  than  the  sum  of  the  number  of  events  (N)  and  the  number 
of  stations  (M).  For  problems  of  interest,  this  commonly  exceeds  200  and  a  systematic 
search  of  a  200  dimensional  solution  space  becomes  very  inefficient.  The  EM  algo¬ 
rithm  is  an  iterative  technique  for  maximizing  the  log-likelihood  function  that  reduces 
to  the  least  squares  solution  in  the  event  that  there  is  no  data  censoring.  A  detailed 
description  of  the  algorithm  is  presented  in  Blandford  and  Shumway  (1982)  and  Bland- 
ford  et  al.  (1984).  The  estimation  procedure  has  been  formalized  in  a  program  called 
MLGLM  (see  Appendix)  written  in  FORTRAN-77  running  on  VAX  11/780,  SUN 
3/160,  Celerity  12600,  and  Cray  XP/II  computer  systems. 

An  additional  problem  with  maximum  likelihood  solutions  is  the  estimation  of 
uncertainties  of  the  estimated  model  parameters.  Although  this  problem  can  be  written 
in  closed  form  for  the  least  squares  solution  of  linear  problems,  the  same  procedure 
would  require  the  determination  of  all  second  order  partial  derivatives  of  the  log- 
likelihood  function.  Calculation  of  such  a  large  Jacobian  for  a  non-linear  function  of 
many  dimensions  is  highly  undesirable.  An  alternative  to  this  is  the  Bootstrap 
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procedure  of  Efron  (1979,  1981).  The  Bootstrapping  procedure  has  been  applied  to  sin¬ 
gle  event  magnitude  estimation  by  McLaughlin  (1988)  for  uncertainty  estimation. 
Using  the  standard  deviations  of  the  perturbing  noise  returned  by  the  EM  algorithm  as 
a  measure  of  the  uncertainty  in  any  single  measurement,  Jih  et  al.  (1988)  proposed 
another  method  to  estimate  the  errors.  In  this  case,  the  o  is  scaled  with  the  number  of 
associated  observations  to  infer  the  precision  of  each  estimated  parameter. 

2.1  Data  Measurement 

The  work  of  Ringdal  (1977)  and  of  von  Seggem  and  Rivers  (1978)  emphasized 
the  importance  in  measuring  or  estimating  noise  and  clipping  levels  in  order  to  be  able 
to  calculate  unbiased  magnitudes.  In  contrast  to  Ringdal  (1986)  or  Lilwall  (1985),  we 
measure  the  non-detection  threshold  levels  and  clipping  levels  for  each  station  and 
event.  It  should,  however,  be  possible  to  extend  this  work  to  estimate  detection  thres¬ 
holds  and/or  clipping  levels  for  use  with  bulletin  data  such  as  the  ISC. 

In  order  to  assure  the  greatest  degree  of  commonality  between  stations  and  events 
in  this  study,  we  have  used  WWSSN  short-period  film  data.  While  use  of  these  data 
has  been  difficult  in  the  past  due  to  their  susceptibility  to  either  clipping  or  masking  by 
system  or  earth  noise,  the  advent  of  the  ability  to  allow  for  noise  and  clipping  within 
the  context  of  a  general  linear  model  reduces  these  difficulties  and  allows  the  WWSSN 
network  to  be  superior  than  others  for  all  except  the  smallest  events.  The  WWSSN 
network  is  valuable  because  it  has  had  a  constant  instrument  response  over  a  long 
period  of  time  with  a  good  spatial  distribution  relative  to  most  test  sites.  The  following 
measurement  procedures  from  Blandford  et  al.  (1984)  are  repeated  here  for  complete- 
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ness. 


2.2  Film  Reading  Procedures  Used  in  this  Study 

1.  Read  the  "a"  (zero-crossing  to  first  peak),  "b"  (first  peak  to  first  trough),  and 
"max"  peak-to-trough  or  trough -to-peak  in  the  first  5  seconds,  in  millimeters. 
Record  readings  in  millimeters  in  the  data  file.  Also,  record  gain  as  written  on 
the  film  and  as  seen  with  particular  viewer. 

2.  Measure  period  for  max  as  peak-to-peak  or  zero-crossing  to  zero-crossing  or 
trough-to-trough  as  in  your  judgement  best  reflects  the  period  of  the  maximum 
energy. 

3.  If  a  weak  signal  must  be  measured  then  try  to  find  a  strong  signal  from  the  same 
test  site  and  by  correlation  try  to  establish  if  a  particular  cycle  is  "a"  or  "b".  In 
general,  it  is  good  practice  to  analyze  events  in  pairs,  with  any  weak  events 
paired  with  a  strong  event. 

4.  If  it  is  not  possible  to  determine  the  exact  location  of  the  "a"  or  "b"  cycles  then 
the  first  clear  down-swing  is  to  be  used  to  place  a  "noise"  limit  on  the  "b"  cycle; 
that  is,  the  true  "b"  amplitude  is  less  than  or  equal  to  that  down-swing.  Apply  a 
similar  method  for  the  "a"  phase  if  there  is  a  clear  up-swing  which  cannot 
confidently  be  said  to  be  a  phase. 

5.  More  commonly,  however,  there  will  be  no  clear  up-swing;  in  this  case  a  noise 
measurement  must  be  made,  and  to  do  so  search  the  preceding  20  seconds  for  the 
largest  peak-to-trough  excursion  in  the  1-2  Hz  frequency  range.  (For  the  a  phase 
find  the  largest  zero- to-peak.)  Often  such  amplitudes  are  less  than  1  mm  on 
WWSSN  film.  The  j)eriod  is  to  be  recorded  as  1  Hz. 

6.  Of  course,  if  no  arrival  at  all  can  be  discerned  then  again  a  noise  measurement 
must  be  made. 

7.  If  the  peak  for  any  phase  can  not  be  found  (clipped)  then  a  conservative  estimate 
of  the  clipping  level  must  be  made.  To  do  this  measure  the  amplitude  of  the 
largest  "turning  point"  visible  and  multiply  by  2.  If  you  can  be  confident  that  the 
largest  turning  point  is  off  the  film  then  you  could  use  the  maximum  distance  to 
the  edge  of  the  film.  If  you  can  not  discern  any  turning  points  and  can  not  deter¬ 
mine  whether  the  trace  is  off  the  film  in  the  first  5  seconds  of  the  P  arrival,  then 
measure  the  largest  on-scale  peak-to-peak  P  coda  that  is  visible.  These  measure¬ 
ments  should  be  recorded  as  clipping  levels. 

8.  In  data  reduction,  all  measurements  with  raw  amplitudes  less  than  1  mm  are  con¬ 
verted  to  noise  or  non-detection  levels.  The  predominant  period,  T,  measured 
from  the  maximum  P  phase  in  the  first  5  seconds,,,,  is  used  to  correct  the  ampli¬ 
tudes  of  "a",  "b",  and  "max"  for  instrument  response,  A(T),  and  to  calculate 
A(T)/T.  For  non-detection  and  clipping,  the  period,  T,  is  assumed  to  be  1  second. 
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3.0  MLE-GLM88 

The  maximum-likelihood  generalized  linear  model,  MLE-GLM88,  for  event  mag¬ 
nitude  estimation  is  performed  with  111  events  as  listed  in  Table  2.  The  mb(a),  and 
mb(b),  and  mb(max)  for  each  of  the  events  are  listed  in  Table  3.  The  "a",  "b",  and 
"max"  magnitudes  for  each  event  are  treated  as  separate  events  in  the  model,  and 
therefore  all  data  from  all  events  for  all  phases  are  used  to  estimate  simultaneously  all 
event  magnitudes  and  station  corrections  in  a  maximum  likelihood  sense. 

The  events  are  grouped  by  test  site:  Cannikin,  Lx>ngshot,  and  Milrow  at  the 
Amchitka  test  site;  Salmon,  Gasbuggy,  Rulison,  Rio  Blanco,  Faultless,  and  Shoal  are 
U.S.  shots  located  away  from  NTS;  Almendro,  Benham,  Bilby,  Bourbon,  Boxcar, 
Cambric,  Chancellor,  Chartreuse,  Corduroy,  Handcar,  Handley,  Harzer,  Kankakee, 
Mast,  Nash,  Piledriver,  Rex,  and  Scotch  at  NTS;  Beryl,  Corundon,  Emeraude,  Grenat, 
Opale,  Rubis,  Saphir,  Tourmaline,  and  Turquoise  at  the  French  southern  Sahara  test 
site;  "dek"  identifies  the  Degelen  Eastern  Kazakhstan  test  site;  "sek"  identifies  the 
Shagan  River  Eastern  Kazakhstan  test  site;  "nnz"  identifies  the  northern  Novaya  Zem- 
lya  test  site;  "snz"  identifies  the  southern  Novaya  Zemlya  test  site;  "pne"  identifies  one 
event  located  in  the  northern  Urals;  "azg"  identifies  events  located  near  Azgir;  "tu" 
identifies  the  Tuomotu  test  site;  "rajl8may74"  is  the  Indian  underground  test  located  in 
Rajasthan,  India;  and  "ch"  identifies  events  located  in  southern  Sinkiang,  China.  The 
number  of  observations  indicates  the  number  of  stations  for  which  a  signal  amplitude, 
noise  level,  or  clipping  level  measurement  was  made.  The  event  magnitude  and  station 
correction  uncertainties  were  estimated  using  the  EM  algorithm.  The  least  squares 
estimates  LS-GLM  are  also  listed  for  comparison.  Although  1 1 1  events  are  listed  for 
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"b"  and  "max"  phase,  5  events  had  no  direct  observations  of  the  "a"  phase  and  only 
106  events  could  be  included  for  the  "a"  phase  magnitude. 
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Table  2.  MLE>GLM88  Event  List 
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Table  2.  (cont’d) 

■rZHn>41Kgfla;Jfcfi»aBrM:ngMBgI«aaiBM[M>ilgft\VS%'^il:<giy^l 


rlrZri?>Mg>>g<!fflTiItMgKl:Trt^Mgl:lEI:MgM:MKV'/ 
■ri:.i  HI  IrJMg  HABiti‘ffglr*:i  Kll  I  a«Mi-f  MW^v/d  I  ;■ 


Teledync  (ieotech 


June  1988 


Final  Report 


Magnitude  Yield  Estimation 


TGAL-87-05 


Table  2.  (cont’d) 
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Table  3a.  MLE-GLM88  mb(a) 


event 

MLE-GLM 

mb 

a 

LS-GLM 

mb 

a 

#  of  observations 
signal  noise  clipped 

cannikin 

6.411 

0.040 

6.416 

0.036 

55 

2 

12 

milrow 

5.947 

0.044 

5.958 

0.037 

51 

5 

0 

longshot 

5.056 

0.038 

5.362 

0.046 

33 

39 

3 

faultless 

5.833 

0.046 

5.841 

0.039 

47 

4 

1 

gasbuggy 

4.144 

0.047 

5.058 

0.188 

2 

47 

0 

no  bianco 

4.042 

0.056 

4.522 

0.153 

3 

32 

0 

rulison 

4.038 

0.048 

4.610 

0.153 

3 

44 

0 

shoal 

4.274 

0.051 

4.847 

0.133 

4 

38 

0 

salmon 

3.450 

0.053 

4.288 

0.266 

1 

38 

0 

almendro 

5.735 

0.062 

5.747 

0.055 

23 

3 

2 

benham 

5.785 

0.046 

5.871 

0.043 

39 

11 

1 

bilby 

5.143 

0.052 

5.380 

0.058 

21 

19 

0 

bourbon 

4.576 

0.047 

4.992 

0.089 

9 

40 

0 

boxcar 

5.833 

0.054 

5.945 

0.052 

26 

10 

1 

cambric 

4.092 

0.049 

4.623 

0.133 

4 

42 

0 

chancellor 

4.882 

0.060 

5.241 

0.100 

7 

22 

1 

chartreuse 

4.872 

0.048 

5.102 

0.064 

17 

30 

1 

chateaugay 

4.467 

0.048 

5.042 

0.153 

3 

42 

2 

corduroy 

4.972 

0.057 

5.180 

0.071 

14 

19 

0 

handcar 

4.305 

0.047 

4.824 

0.119 

5 

45 

0 

handley 

6.064 

0.050 

6.054 

0.041 

42 

1 

0 

harzer 

5.005 

0.053 

5.362 

0.077 

12 

25 

1 

kankakee 

4.334 

0.046 

4.976 

0.108 

6 

45 

0 

mast 

5.406 

0.059 

5.527 

0.054 

24 

7 

0 

nash 

4.750 

0.045 

5.119 

0.069 

15 

38 

0 

piledriver 

4.918 

0.046 

5.094 

0.055 

23 

29 

0 

rex 

3.867 

0.045 

4.451 

0.188 

2 

51 

0 

scotch 

5.068 

0.048 

5.215 

0.058 

21 

26 

1 

azgl4iul79 

4.841 

0.099 

5.235 

0.153 

3 

7 

1 

azgl7ian79 

5.528 

0.088 

5.602 

0.089 

9 

2 

3 

azgl7oct78 

5.287 

0.095 

5.494 

0.089 

9 

2 

1 

azgl8dec78 

5.379 

0.095 

5.682 

0.108 

6 

4 

2 

azg22apr66 

3.863 

0.091 

5.009 

0.266 

1 

12 

0 

azg22dec7 I 

5.483 

0.088 

5.629 

0.074 

13 

1 

0 

azg24oct79 

4.863 

0.110 

5.251 

0.188 

2 

6 

1 

azg29iul76 

5.105 

0.045 

5.444 

0.055 

23 

26 

4 

azg30sep77 

4.043 

0.045 

5.025 

0.188 

2 

51 

0 
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TGAL-87.05 

event 

MLE-GLM 

mb 

a 

LS-GLM 

mb 

a 

#  of  observations 
signal  noise  clipped 

dekllsep69 

3.976 

0.043 

4.761 

0.153 

3 

57 

0 

dekl8dec66 

5.253 

0.040 

5.459 

0.041 

42 

25 

0 

dek22apr78 

4.477 

0.060 

4.892 

0.094 

8 

22 

0 

dek25apr71 

5.278 

0.050 

5.459 

0.049 

30 

13 

0 

dek26feb67 

5.369 

0.041 

5.510 

0.038 

48 

18 

0 

dek26mar78 

4.965 

0.058 

5.235 

0.069 

15 

17 

0 

dek29mar77 

4.315 

0.054 

4.771 

0.094 

8 

29 

0 

dek30dec71 

5.018 

0.074 

5.263 

0.077 

12 

8 

0 

dek30jul77 

4.196 

0.053 

4.816 

0.133 

4 

34 

0 

nnz01oct81 

5.230 

0.046 

5.360 

0.045 

35 

14 

3 

nnz01sep77 

5.093 

0.060 

5.270 

0.057 

22 

8 

0 

nnz07nov68 

5.597 

0.041 

5.695 

0.037 

51 

11 

4 

nnzl0aug78 

5.395 

0.043 

5.447 

0.040 

45 

8 

7 

nnzlloctSO 

5.185 

0.046 

5.256 

0.046 

33 

14 

5 

nnzl2sep73 

6.369 

0.049 

6.325 

0.050 

28 

0 

17 

nnzl4oct69 

5.762 

0.040 

5.809 

0.037 

53 

9 

6 

nnzi4oct70 

6.434 

0.043 

6.424 

0.039 

47 

1 

10 

nnzl8aug83 

5.327 

0.052 

5.414 

0.049 

30 

8 

2 

nnz20oct76 

4.027 

0.043 

4.641 

0.108 

6 

54 

0 

nnz21oct67 

5.397 

0.042 

5.538 

0.039 

46 

15 

1 

nnz21oct75 

6.106 

0.051 

6.112 

0.048 

31 

1 

9 

nnz23aug75 

6.119 

0.051 

6.134 

0.045 

35 

0 

6 

nnz25oct84 

5.158 

0.061 

5.318 

0.063 

18 

10 

1 

nnz27oct66 

6.071 

0.039 

6.061 

0.035 

58 

2 

11 

nnz27sep71 

6.269 

0.049 

6.248 

0.048 

31 

1 

13 

nnz28aug72 

5.993 

0.050 

5.996 

0.046 

34 

2 

8 

nnz29aug74 

6.128 

0.050 

6.155 

0.047 

32 

1 

10 

snz02nov74 

6.504 

0.051 

6.440 

0.053 

25 

0 

16 

snzl8oct75 

6.227 

0.051 

6.228 

0.047 

32 

1 

9 

snz27oc73a 

6.650 

0.053 

6.555 

0.063 

18 

0 

20 

snz27oc73c 

3.522 

0.053 

4.216 

0.266 

1 

38 

0 

snz27sep73 

5.182 

0.054 

5.326 

0.049 

29 

8 

0 

sek04juI76 

5.245 

0.074 

5.422 

0.069 

15 

5 

0 

sek07dec76 

4.990 

0.076 

5.451 

0.084 

10 

9 

0 

sekl  ljun78 

5.272 

0.076 

5.452 

0.071 

14 

5 

0 

sekl4dec73 

5.242 

0.042 

5.476 

0.048 

31 

29 

1 

sekl4sep80 

5.436 

0.056 

5.654 

0.053 

25 

10 

0 
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Magnitude  Yield  Estimation 


TGAL-87-05 


Table  3b.  MLE-GLM88  1115(b) 


event 

MLE-GLM 

nib 

0 

LS-GLM 

nib 

0 

#  of  observations 
signal  noise  clipped 

cannikin 

6.667 

0.040 

6.598 

0.037 

51 

0 

18 

longshot 

5.432 

0.038 

5.537 

0.036 

54 

18 

3 

milrow 

6.198 

0.044 

6.146 

0.036 

54 

0 

2 

faultless 

6.158 

0.046 

6.128 

0.038 

48 

1 

3 

gasbuggy 

4.404 

0.047 

4.890 

0.089 

9 

40 

0 

no  bianco 

4.537 

0.056 

4.903 

0.084 

10 

25 

0 

nilison 

4.151 

0.048 

4.676 

0.119 

5 

42 

0 

salmon 

3.979 

0.053 

4.387 

0.119 

5 

34 

0 

shoal 

4.426 

0.051 

4.775 

0.084 

10 

32 

0 

almendro 

6.031 

0.062 

6.005 

0.052 

26 

0 

2 

benham 

6.103 

0.046 

6.056 

0.039 

46 

1 

4 

bilby 

5.407 

0.052 

5.459 

0.046 

34 

6 

0 

bourbon 

4.711 

0.047 

5.014 

0.071 

14 

35 

0 

boxcar 

6.178 

0.054 

6.108 

0.046 

33 

0 

4 

cambric 

4.341 

0.049 

4.660 

0.084 

10 

36 

0 

chancellor 

5.176 

0.060 

5.313 

0.066 

16 

13 

1 

chartreuse 

5.001 

0.048 

5.135 

0.053 

25 

22 

1 

chateaugay 

4.882 

0.048 

5.198 

0.071 

14 

30 

3 

corduroy 

5.091 

0.057 

5.359 

0.069 

15 

18 

0 

handcar 

4.496 

0.047 

4.719 

0.071 

14 

36 

0 

handley 

6.309 

0.050 

6.298 

0.041 

41 

1 

1 

harzer 

5.305 

0.053 

5.388 

0.051 

27 

10 

1 

kankakee 

4.591 

0.046 

4.837 

0.061 

19 

32 

0 

mast 

5.735 

0.059 

5.761 

0.049 

29 

2 

0 

nash 

4.914 

0.045 

5.085 

0.052 

26 

27 

0 

piledriver 

5.191 

0.046 

5.283 

0.044 

36 

16 

0 

rex 

4.372 

0.045 

4.669 

0.077 

12 

41 

0 

scotch 

5.348 

0.048 

5.470 

0.047 

32 

15 

1 

azgl4jul79 

5.375 

0.099 

5.509 

0.094 

8 

1 

2 

azgl7jan79 

5.880 

0.088 

5.885 

0.084 

10 

0 

4 

azgl7oct78 

5.737 

0.095 

5.781 

0.094 

8 

0 

4 

azgl8dec78 

5.751 

0.095 

5.859 

0.089 

9 

1 

2 

azg22apr66 

4.100 

0.091 

4.494 

0.153 

3 

10 

0 

azg22dec7 1 

5.835 

0.088 

5.879 

0.074 

13 

0 

1 

azg24oct79 

5.696 

0.110 

5.691 

0.153 

3 

1 

5 

azg25apr75 

3.915 

0.080 

4.817 

0.266 

1 

16 

0 

azg29jul76 

5.580 

0.045 

5.612 

0.042 

40 

6 

7 

azg30sep77 

4.579 

0.045 

4.976 

0.069 

15 

37 

1 
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Table  3b.  (cont’d) 


event 

MLE-GLM 

nib 

a 

LS-GLM 

nib 

o 

#  of  observations 
signal  noise  clipped 

dekllsep69 

4.232 

0.043 

4.710 

0.089 

9 

51 

0 

dekl8dec66 

5.491 

0.040 

5.608 

0.037 

53 

14 

0 

dek22apr78 

4.781 

0.060 

5.051 

0.064 

17 

13 

0 

dek25apr71 

5.552 

0.050 

5.629 

0.043 

38 

5 

0 

dek26feb67 

5.609 

0.041 

5.705 

0.037 

52 

12 

2 

dek26mar78 

5.284 

0.058 

5.392 

0.055 

23 

9 

0 

dek29mar77 

4.710 

0.054 

4.953 

0.061 

19 

18 

0 

dek30dec71 

5.366 

0.074 

5.452 

0.064 

17 

3 

0 

dek30jul77 

4.595 

0.053 

4.918 

0.071 

14 

24 

0 

nnzOloctSl 

5.492 

0.046 

5.511 

0.041 

42 

4 

6 

nnz01sep77 

5.426 

0.060 

5.492 

0.052 

26 

3 

1 

nnz07nov68 

5.845 

0.041 

5.842 

0.035 

57 

2 

7 

nnzl0aug78 

5.630 

0.043 

5.590 

0.041 

41 

3 

16 

nnzlloct80 

5.445 

0.046 

5.472 

0.041 

42 

5 

5 

nnzl2sep73 

6.699 

0.049 

6.615 

0.058 

21 

0 

24 

nnzl4oct69 

5.960 

0.040 

5.979 

0.036 

55 

5 

8 

nnzl4oct70 

6.648 

0.043 

6.546 

0.043 

39 

0 

19 

nnzl8aug83 

5.531 

0.052 

5.567 

0.049 

29 

6 

5 

nnz20oct76 

4.353 

0.043 

4.717 

0.066 

16 

43 

1 

nnz21oct67 

5.595 

0.042 

5.685 

0.037 

52 

9 

1 

nnz21oct75 

6.348 

0.051 

6.302 

0.049 

29 

0 

12 

nnz23aug75 

6.381 

0.051 

6.339 

0.049 

30 

0 

11 

nnz25oct84 

5.434 

0.061 

5.452 

0.055 

23 

3 

3 

nnz27oct66 

6.306 

0.039 

6.255 

0.036 

54 

1 

16 

nnz27sep71 

6.487 

0.049 

6.401 

0.052 

26 

0 

19 

nnz28aug72 

6.254 

0.050 

6.200 

0.047 

32 

0 

12 

nnz29aug74 

6.397 

0.050 

6.312 

0.053 

25 

1 

17 

snz02nov74 

6.799 

0.051 

6.633 

0.064 

17 

0 

24 

snzl8oct75 

6.519 

0.051 

6.427 

0.052 

26 

0 

16 

snz27oc73a 

6.879 

0.053 

6.698 

0.071 

14 

0 

24 

snz27oc73b 

4.007 

0.054 

4.340 

0.108 

6 

31 

0 

snz27oc73c 

3.854 

0.053 

4.335 

0.133 

4 

35 

0 

snz27sep73 

5.474 

0.054 

5.570 

0.048 

31 

5 

1 

sek04jul76 

5.609 

0.074 

5.690 

0.063 

18 

1 

1 

sek07dec76 

5.439 

0.078 

5.485 

0.063 

18 

0 

0 

sekl  ljun78 

5.579 

0.078 

5.638 

0.064 

17 

1 

0 

sekl4dec73 

5.546 

0.042 

5.659 

0.041 

41 

16 

4 
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Table  3b.  (cont’d) 

MLE-GLM  LS-GLM  #  of  observations 


event 

mb 

a 

mb 

a 

signal 

noise 

clipped 

sekMsepSO 

5.748 

0.056 

5.861 

0.049 

30 

5 

0 

sekl5jan65 

5.736 

0.046 

5.763 

0.038 

48 

1 

3 

sekl5sep78 

5.708 

0.054 

5.693 

0.048 

31 

2 

4 

sekl9jun68 

4.999 

0.056 

5.104 

0.052 

26 

7 

1 

sek23juI73 

5.992 

0.051 

6.041 

0.042 

40 

1 

1 

sek23jun79 

5.852 

0.050 

5.914 

0.043 

39 

4 

0 

sek27apr75 

5.288 

0.076 

5.423 

0.064 

17 

2 

0 

sek3(}nov69 

5.757 

0.056 

5.803 

0.046 

33 

1 

0 

pne29aug74 

4.401 

0.049 

4.637 

0.059 

20 

25 

0 

ch03oct84 

4.756 

0.070 

5.016 

0.094 

8 

14 

0 

ch06oct83 

5.030 

0.059 

5.232 

0.071 

14 

16 

1 

chl7oct76 

4.142 

0.048 

4.745 

0.119 

5 

42 

0 

chl9dec84 

3.954 

0.091 

4.663 

0.266 

1 

12 

0 

ch22sep69 

4.747 

0.051 

4.866 

0.053 

25 

17 

0 

ch27oct75 

4.402 

0.054 

4.672 

0.084 

10 

27 

0 

rajl8may74 

4.307 

0.060 

4.746 

0.119 

5 

25 

0 

tu03dec80 

4.982 

0.050 

5.197 

0.054 

24 

19 

0 

tul9apr83 

5.193 

0.069 

5.304 

0.063 

18 

5 

0 

tul9feb77 

4.363 

0.049 

4.700 

0.080 

11 

34 

0 

tul9jul80 

5.160 

0.051 

5.279 

0.047 

32 

9 

1 

tul9mar77 

5.445 

0.063 

5.534 

0.058 

21 

6 

0 

tu23mar80 

5.106 

0.050 

5.243 

0.057 

22 

19 

3 

tu24nov77 

5.364 

0.060 

5.488 

0.053 

25 

5 

0 

tu25jul79 

5.563 

0.078 

5.632 

0.066 

16 

2 

0 

tu25jul82 

5.035 

0.056 

5.200 

0.058 

21 

14 

0 

tu25may83 

5.446 

0.078 

5.458 

0.069 

15 

3 

0 

tu30nov78 

5.235 

0.047 

5.384 

0.048 

31 

17 

1 

beryl 

4.755 

0.078 

5.028 

0.084 

10 

8 

0 

corundon 

3.897 

0.045 

4.500 

0.133 

4 

50 

0 

emeraude 

4.218 

0.051 

4.795 

0.100 

7 

34 

0 

grenat 

4.495 

0.041 

4.845 

0.057 

22 

43 

0 

opale 

3.847 

0.045 

4.507 

0.153 

3 

50 

0 

rubis 

5.155 

0.048 

5.316 

0.044 

36 

12 

0 

saphir 

5.460 

0.041 

5.563 

0.036 

54 

8 

1 

tourmaline 

4.424 

0.040 

4.821 

0.063 

18 

49 

0 

turqoise 

3.925 

0.040 

4.492 

0.133 

4 

63 

0 
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Table  3c.  MLE-GLM88  in|,(max) 


event 

MLE-GLM 

mb 

o 

LS-GLM 

mb 

a 

#  of  observations 
signal  noise  clipped 

cannikin 

6.916 

0.040 

6.832 

0.038 

49 

0 

20 

milrow 

6.497 

0.044 

6.433 

0.037 

52 

0 

4 

longshot 

5.823 

0.038 

5.854 

0.032 

67 

4 

3 

faultless 

6.458 

0.046 

6.433 

0.038 

48 

1 

3 

gasbuggy 

4.654 

0.047 

5.177 

0.080 

11 

38 

0 

no  bianco 

4.803 

0.056 

5.042 

0.069 

15 

20 

0 

rulison 

4.524 

0.048 

4.915 

0.084 

10 

37 

0 

shoal 

4.713 

0.051 

5.047 

0.071 

14 

28 

0 

salmon 

4.183 

0.053 

4.503 

0.108 

6 

33 

0 

almendro 

6.238 

0.062 

6.215 

0.052 

26 

0 

2 

benham 

6.357 

0.046 

6.285 

0.041 

43 

1 

7 

bilby 

5.655 

0.052 

5.666 

0.044 

37 

3 

0 

bourbon 

4.897 

0.047 

5.199 

0.063 

18 

31 

0 

boxcar 

6.399 

0.054 

6.329 

0.046 

33 

0 

4 

cambric 

4.550 

0.049 

4.901 

0.077 

12 

34 

0 

chancellor 

5.329 

0.060 

5.503 

0.066 

16 

13 

1 

chartreuse 

5.242 

0.048 

5.338 

0.048 

31 

16 

1 

chateaugay 

5.063 

0.048 

5.399 

0.064 

17 

28 

2 

corduroy 

5.282 

0.057 

5.527 

0.063 

18 

15 

0 

handcar 

4.629 

0.047 

4.845 

0.066 

16 

34 

0 

handley 

6.481 

0.050 

6.475 

0.041 

41 

1 

1 

harzer 

5.527 

0.053 

5.545 

0.048 

31 

6 

1 

kankakee 

4.842 

0.046 

5.045 

0.054 

24 

27 

0 

mast 

5.971 

0.059 

5.990 

0.049 

30 

1 

0 

nash 

5.144 

0.045 

5.285 

0.048 

31 

22 

0 

piledriver 

5.431 

0.046 

5.500 

0.043 

38 

13 

1 

rex 

4.717 

0.045 

4.933 

0.066 

16 

36 

1 

scotch 

5.603 

0.048 

5.655 

0.043 

39 

8 

1 

azgl4jul79 

5.705 

0.099 

5.747 

0.084 

10 

0 

1 

azgl7jan79 

6.161 

0.088 

6.192 

0.084 

10 

0 

4 

azgl7oct78 

6.112 

0.095 

6.122 

0.100 

7 

0 

5 

azgl8dec78 

6.118 

0.095 

6.161 

0.089 

9 

0 

3 

azg22apr66 

4.181 

0.091 

4.657 

0.153 

3 

10 

0 

azg22dec71 

6.173 

0.088 

6.208 

0.077 

12 

0 

2 

azg24oct79 

5.980 

O.llO 

5.961 

0.153 

3 

0 

6 

azg25apr75 

3.954 

0.080 

4.857 

0.266 

1 

16 

0 

azg29jul76 

5.865 

0.045 

5.891 

0.041 

41 

5 

7 

azg30sep77 

4.816 

0.045 

5.141 

0.058 

21 

31 

1 
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Table  3c.  (cont’d) 

MLE-GLM  LS-GLM  #  of  observations 


event 

mb 

a 

mb 

a 

signal 

noise 

clipped 

dekllsep69 

4.571 

0.043 

4.932 

0.061 

19 

41 

0 

dek;18dec66 

5.704 

0.040 

5.764 

0.035 

58 

8 

1 

dek22apr78 

5.028 

0.060 

5.196 

0.058 

21 

9 

0 

dek25apr71 

5.739 

0.050 

5.831 

0.043 

38 

5 

0 

dek26feb67 

5.829 

0.041 

5.892 

0.037 

51 

9 

6 

dek26mar78 

5.505 

0.058 

5.581 

0.052 

26 

6 

0 

dek29mar77 

4.989 

0.053 

5.172 

0.053 

25 

14 

0 

dek30dec71 

5.544 

0.074 

5.639 

0.064 

17 

3 

0 

dek30jul77 

4.859 

0.053 

5.069 

0.057 

22 

16 

0 

nnz01oct81 

5.651 

0.046 

5.678 

0.041 

43 

4 

5 

nnz01sep77 

5.574 

0.059 

5.610 

0.051 

27 

2 

2 

nnz07nov68 

6.040 

0.041 

6.043 

0.034 

60 

1 

5 

nnzl0aug78 

5.861 

0.043 

5.812 

0.043 

39 

3 

18 

nnzlloct80 

5.661 

0.046 

5.685 

0.041 

42 

4 

6 

nnzl2sep73 

6.795 

0.049 

6.715 

0.054 

24 

0 

21 

nnzl4oct69 

6.131 

0.040 

6.117 

0.035 

59 

2 

7 

nnzl4oct70 

6.831 

0.043 

6.709 

0.044 

36 

0 

22 

nnzl8aug83 

5.709 

0.052 

5.738 

0.049 

30 

5 

5 

nnz20oct76 

4.660 

0.043 

4.913 

0.053 

25 

34 

1 

nnz21oct67 

5.771 

0.042 

5.787 

0.036 

54 

5 

3 

nnz21oct75 

6.561 

0.051 

6.445 

0.054 

24 

0 

17 

nnz23aug75 

6.506 

0.051 

6.458 

0.049 

29 

0 

12 

nnz25oct84 

5.606 

0.061 

5.608 

0.057 

22 

3 

4 

nnz27oct66 

6.449 

0.039 

6.397 

0.035 

57 

0 

14 

nnz27sep71 

6.635 

0.049 

6.515 

0.054 

24 

0 

21 

nnz28aug72 

6.382 

0.050 

6.316 

0.046 

33 

0 

11 

nnz29aug74 

6.582 

0.050 

6.495 

0.053 

25 

0 

18 

snz02nov74 

7.021 

0.051 

6.792 

0.077 

12 

0 

29 

snzl8oct75 

6.836 

0.051 

6.716 

0.058 

21 

0 

21 

snz27oc73a 

7.099 

0.053 

6.961 

0.071 

14 

0 

24 

snz27oc73b 

4.156 

0.054 

4.454 

0.089 

9 

28 

0 

snz27oc73c 

3.875 

0.053 

4.383 

0.133 

4 

35 

0 

snz27sep73 

5.717 

0.054 

5.783 

0.046 

33 

3 

1 

sek04jul76 

5.928 

0.074 

5.904 

0.069 

15 

0 

5 

sek07dec76 

5.614 

0.072 

5.760 

0.063 

18 

2 

1 

seklljun78 

5.878 

0.076 

5.880 

0.063 

18 

0 

1 

sekl4dec73 

5.773 

0.042 

5.819 

0.039 

47 

8 

6 
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_ Table  3c.  (cont’d) _ 

MLE-GLM  LS-GLM  #  of  observations 


event 

mb 

a 

ntjj 

a 

signal 

noise 

clipped 

sekl4sep80 

5.985 

0.052 

6.092 

0.049 

29 

5 

6 

sekl5jan65 

5.883 

0.046 

5.917 

0.038 

49 

1 

2 

sekl5sep78 

5.899 

0.054 

5.866 

0.048 

31 

1 

5 

sekl9jun68 

5.255 

0.056 

5.308 

0.049 

29 

3 

2 

sek23jul73 

6.177 

0.051 

6.228 

0.042 

40 

1 

1 

sek23jun79 

6.050 

0.049 

6.101 

0.042 

40 

3 

3 

sek27apr75 

5.531 

0.072 

5.614 

0.061 

19 

1 

1 

sek30nov69 

5.958 

0.056 

5.987 

0.046 

34 

0 

0 

pne29aug74 

4.726 

0.049 

4.889 

0.051 

27 

18 

0 

ch03oct84 

5.005 

0.070 

5.227 

0.084 

10 

12 

0 

ch06oct83 

5.245 

0.059 

5.393 

0.064 

17 

13 

1 

chl7oct76 

4.532 

0.048 

4.848 

0.074 

13 

34 

0 

chl9dec84 

4.366 

0.091 

4.825 

0.153 

3 

10 

0 

ch22sep69 

5.134 

0.051 

5.183 

0.049 

30 

12 

0 

ch27oct75 

4.590 

0.054 

4.877 

0.077 

12 

25 

0 

rajl8may74 

4.565 

0.060 

5.047 

0.100 

7 

23 

0 

tu03dec80 

5.333 

0.050 

5.449 

0.047 

32 

11 

0 

tul9apr83 

5.491 

0.069 

5.495 

0.057 

22 

1 

0 

tul9feb77 

4.608 

0.049 

4.924 

0.066 

16 

29 

0 

tul9juI80 

5.515 

0.051 

5.510 

0.043 

38 

2 

2 

tul9mar77 

5.646 

0.063 

5.759 

0.059 

20 

6 

1 

tu23mar80 

5.358 

0.050 

5.470 

0.051 

27 

14 

3 

tu24nov77 

5.659 

0.057 

5.652 

0.046 

33 

0 

0 

tu25jul79 

5.856 

0.078 

5.838 

0.063 

18 

0 

0 

tu25jul82 

5.211 

0.056 

5.406 

0.057 

22 

13 

0 

tu25may83 

5.778 

0.078 

5.721 

0.063 

18 

0 

0 

tu30nov78 

5.611 

0.047 

5.668 

0.042 

40 

7 

2 

beryl 

4.995 

0.078 

5.237 

0.077 

12 

6 

0 

corundon 

4.212 

0.045 

4.661 

0.080 

11 

43 

0 

emeraude 

4.541 

0.052 

4.830 

0.069 

15 

25 

0 

grenat 

4.766 

0.041 

5.025 

0.047 

32 

32 

1 

opale 

3.889 

0.045 

4.603 

0.153 

3 

50 

0 

rubis 

5.433 

0.048 

5.517 

0.041 

42 

5 

1 

saphir 

5.729 

0.041 

5.769 

0.037 

53 

5 

5 

tourmaline 

4.643 

0.040 

4.976 

0.051 

27 

40 

0 

turqoise 

4.227 

0.040 

4.689 

0.077 

12 

55 

0 
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3.1  Censoring  versus  No-censoring 

Figure  1  plots  the  GLM  versus  the  MLE-GLM  magnitudes  for  the  complete  set  of 
events  for  the  "a",  "b",  and  "max"  phase.  The  results  show  that  in  using  the  stations 
belonging  to  the  WWSSN,  there  is  only  a  small  magnitude  range  from  about  5.7  to  6.5 
where  ignoring  the  censoring  effects  of  clipping  or  non-detection  does  not  lead  to 
biased  magnitudes.  For  events  less  than  mj,  of  5,  the  GLM  and  MLE-GLM  magnitude 
bias  for  the  "a"  phase  may  be  up  to  a  full  magnitude  unit.  For  the  "b"  and  "max" 
phases,  the  magnitude  bias  between  the  GLM  and  MLE-GLM  averages  to  be  about 
0.6. 

If  all  calibration  events  could  be  selected  in  the  magnitude  range  of  5.7  to  6.5 
where  no  bias  is  found  between  GLM  and  MLE-GLM  and  if  those  calibration  events 
were  truly  representative  of  the  test  site,  then  an  MLE-GLM  would  not  be  necessary. 
However,  in  order  to  have  a  representative  sample  of  explosions  from  a  test  site  and  a 
representative  sample  from  all  test  sites  for  comparison,  it  rapidly  becomes  necessary 
to  select  events  from  outside  this  range  of  magnitudes.  An  example  is  the  three  events 
at  Amchitka,  where  the  WWSSN  stations  underwent  changes  in  network  coverage  in 
the  late  1960’s  between  the  detonations  of  Longshot,  Milrow,  and  Cannikin.  This 
reconfiguration  of  the  network  coverage  may  produce  inconsistency  in  the  m),  estima¬ 
tion  which  results  in  Longshot  not  being  an  adequate  calibration  event  for  Milrow  and 
Cannikin. 
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mb(a) 


m^(a)  (no  — censoring) 


Figure  la.  nib(a)  event  magnitudes  from  LS-GLM  (no-censoring)  estimates  plotted 
against  MLE-GLM  (censoring)  estimates  for  106  explosions.  Note  that  there  is  only  a 
narrow  range  of  magnitudes  (5.7  to  6.5)  for  which  the  two  sets  of  magnitudes  are  on 
average  the  same.  The  LS-GLM  event  magnitudes  are  too  large  at  the  lower  m^’s. 
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Figure  lb.  mj,(b)  event  magnitudes  from  LS-GLM  (no-censoring)  estimates  plotted 
against  MLE-GLM  (censoring)  estimates  for  109  explosions.  Note  that  there  is  only  a 
narrow  range  of  magnitudes  (5.5  to  6.5)  for  which  the  two  sets  of  magnitudes  are  on 
average  the  same.  The  LS-GLM  event  magnitudes  are  too  large  at  the  lower  m^’s  and 
too  small  at  the  higher  m^’s. 
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m,  (max) 


6.5 


cn 

.5  6.0 

!>_ 

o 

CO 

£  5.5 

u 


^  5.0 

o 

£ 


4.5 


4.0 


5.5 


4.0  5.0  6.0  7.0 

m^(max)  (no  — censoring) 

Figure  Ic.  mb(max)  event  magnitudes  from  LS-GLM  (no-censoring)  estimates  plotted 
against  MLE-GLM  (censoring)  estimates  for  109  explosions.  Note  that  there  is  only  a 
narrow  range  of  magnitudes  (5.5  to  6.5)  for  which  the  two  sets  of  magnitudes  are  on 
average  the  same.  The  LS-GLM  event  magnitudes  are  too  large  at  the  lower  m^’s  and 
too  small  at  the  higher  m^’s. 
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3.2  Phase  Comparison 

The  GLM  and  MLE-GLM  mj,  magnitudes  estimated  for  all  events  using  both  the 
censoring  and  non-censoring  information  for  the  three  phases  totaling  32,328  event- 
phases  readings  are  plotted  in  Figures  2  and  3  respectively.  In  the  case  of  GLM  with 
no-censoring  information,  the  mt,(a)  values  are  on  average  0.4  magnitude  unit  less  than 
m(,(max)  and  0.16  magnitude  unit  less  than  mb(b)  for  mb(a)  over  5.3.  The  m^  bias 
between  the  m^Cb)  and  mb(max)  is  0.2.  For  events  less  than  5.3,  the  m,,  bias 
decreases.  The  departure  from  the  normal  magnitude  bias  for  the  mb(a)  m^Cb)  and 
m(,(a)  vs  m(,(max)  may  be  due  to  the  low  signal  detection  levels  for  some  of  the  lower 
yield  events. 

In  the  case  of  MLE-GLM  when  censoring  data  are  included,  the  mj,  biases 
between  mb(a)  and  other  phases  have  increased  compared  to  those  without  including 
the  censoring  information  (Figure  3).  The  mb  bias  is  approximately  0.28  for  mb(a)  vs 
mb(b),  0.19  for  mb(b)  vs  mb(max),  and  0.48  for  mb(a)  vs  mb(max).  The  inclusion  of 
censoring  information  provides  better  constraints  on  the  mb’s  of  the  low  yield  events 
and  therefore  furnishes  a  more  coherent  set  of  magnitude  relationship  among  the 
phases. 

3.3  Multiple  versus  Single  Phase  Determination 

The  maximum-likelihood  mb  estimations  arc  performed  in  two  different  data 
configurations.  The  first  one  is  to  perform  LS-GLM  and  MLE-GLM  estimates  using 
data  for  all  three  phases  for  111  events  assuming  they  are  correlated.  The  other  esti¬ 
mates  are  obtained  by  grouping  the  data  for  each  individual  phase  for  all  events  and 
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m^(no- censoring) 


Figure  2b.  LS-GUM  (no-censoring)  m(,’s  estimates  for  the  "a"  phase  plotted  against  the 
"max"  phase.  The  mb(max)’s  are  on  the  average  0.4  magnitude  unit  greater  than  the 
mb(a)’s  for  the  larger  events. 
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6.5 


m,(  no  — censoring) 


6.0 


X  5  5 


5.0 


4.5 


4.0 


3.5 


4.0 


5.0 


mh(b) 


Figure  2c.  LS-GLM  (no-censoring)  m|,’s  estimates  for  the  "b"  phase  plotted  against  the 
"max"  phase.  The  mb(max)’s  are  on  the  average  0.2  magnitude  unit  greater  than  the 
mb(a)’s  for  the  larger  events. 
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m,  (censoring) 


6.0 


5.5 


5.0 


4.5 


4.0  h  . 


3.5 


4.0 


5.0 

mbCa) 


7.0 


Figure  3a.  MLE-GLM  mb’s  estimates  for  the  "a"  phase  plotted  against  those  for  the 
"b"  phase.  These  mb’s  are  derived  using  all  three  phases.  The  mb(b)’s  are  on  the  aver¬ 
age  0.28  magnitude  unit  greater  than  the  mb(a)’s. 


Teledyne  Geotech 


June  1988 


Final  report 


Magnitude  Yield  Estimation 


TGAL-87-05 
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5.0 


Figure  3c.  MLE-GLM  m^’s  estimates  for  the  "b"  phase  plotted  against  those  for  the 
"max"  phase.  These  m^’s  are  derived  using  all  three  phases.  The  mb(max)’s  are  on  the 
average  0.19  magnitude  unit  greater  than  the  mb(a)’s. 
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performing  three  separate  inversions.  The  results  are  plotted  in  Figures  4  and  5  for  the 
LS-GLM  and  MLE-GLM  respectively.  When  no  censoring  information  is  taken  into 
account,  the  mb(a)  estimates  using  just  the  "a"  phases  alone  are  on  the  average  higher 
than  the  mb(a)’s  obtained  by  inverting  all  three  phases  together  (Figure  4a).  On  the 
other  hand,  the  mi,(b)’s  obtained  using  all  phases  appear  to  be  slightly  lower  than  those 
obtained  using  just  the  "b"  phase  (Figure  4b).  There  is  no  apparent  bias  between  the 
mb(max)’s  obtained  using  all  phases  and  just  the  "max"  phase. 

When  censored  data  are  taken  into  account,  the  mb(a)’s  obtained  using  all  phases 
appears  to  be  slightly  lower  than  those  obtained  by  using  just  the  "a"  phase  for  the 
lower  yield  event  (Figure  5a).  In  the  case  of  mb(max),  the  magnitude  estimates  using 
just  the  "max"  phase  are  slightly  higher  than  those  obtained  using  all  three  phases  for 
the  lower  yield  events  (Figure  5c),  inverse  to  the  pattern  for  mb(a)’s.  No  apparent  bias 
is  observed  between  the  two  sets  of  mb(b)’s.  But  in  all  cases,  the  mb  biases  for  esti¬ 
mates  using  all  phases  versus  those  using  the  individual  phases  are  small,  being  no 
more  than  0.03  magnitude  unit. 

3.4  Station  Effects 

The  estimation  of  the  station  terms  has  a  controlling  effect  on  the  estimation  of 
mb’s.  The  station  effect  is  a  function  of  the  local  elastic  and  anelastic  structure  near 
the  station.  A  secondary  dependence  may  be  attributed  to  the  azimuthal  and  incidence 
angle  of  the  ray  paths,  due  to  a  focusing-defocusing  effect  near  the  station.  Stations 
corrections  have  been  constructed  by  several  researchers  using  various  sets  of  data. 
Ringdal  (1984)  derived  station  corrections  using  nearly  40,000  P-wave  amplitudes. 
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no  — censoring 

m(^(nnQx)  (all  phases) 

Figure  4c.  LS-GLM  mj,(max)  estimates  derived  from  using  just  the  "max"  phase  plot¬ 
ted  versus  those  derived  from  using  all  three  phases.  These  two  sets  of  magnitudes 
agree  very  well. 
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m^(a)  (all  phases) 

Figure  5a.  MLE-GLM  mj,(a)  estimates  derived  from  using  just  the  "a"  phase  plotted 
versus  those  derived  from  using  all  three  phases.  There  is  a  slight  indication  that  the 
m(,(a)’s  for  just  the  "a"  phase  zue  biased  low  towards  low  m^’s. 
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Figure  5b.  MLE-GLM  m|,(a)  estimates  derived  from  using  just  the  "b"  phase  plotted 
versus  those  derived  from  using  all  three  phases.  These  two  sets  of  magnitudes  agree 
very  well. 
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Figure  5c.  MLE-GLM  mb(max)  estimates  derived  from  using  just  the  "max"  phase 
plotted  versus  those  derived  from  using  all  three  phases.  These  two  sets  of  magnitudes 
agree  very  well. 
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whereas  Blandford  et  al.  (1984)  used  data  from  nuclear  test  sites  to  derive  average  sta¬ 
tion  corrections. 

Table  4  lists  the  maximum  likelihood  station  effect  estimates  and  their  uncertain¬ 
ties,  a,  estimated  from  a  bootstrap  procedure.  The  number  of  observations,  noise  level 
estimates,  and  clipping  level  estimates  at  each  station  are  also  listed.  The  RMS  station 
standard  deviation  is  0.0279,  and  the  RMS  station  effect  is  0.2112.  We  have  assumed 
that  the  station  corrections  for  "a",  "b",  and  "max"  phases  are  the  same.  Due  to  the 
fact  that  the  station  effects  were  determined  in  a  maximum  likelihood  sense,  the  cen¬ 
soring  information,  such  as  non-detection  or  clipping,  levels  was  also  taken  into 
account.  Stations  that  may  have  often  clipped  or  not  detected  signals  were  assigned 
appropriate  station  corrections  that  reflected  their  tendency  to  have  larger  or  smaller 
amplitudes. 

The  station  effects  pertaining  to  e.  ,n  phase  are  studied  separately  by  not  making 
the  assumption  that  the  station  effects  are  the  same  for  the  different  phases.  Figure  6 
plots  the  MLE-GLM  station  effects  derived  from  each  of  the  three  phases  separately. 
The  station  effects  range  from  -0.8  to  0.6  and  show  quite  a  scatter  among  the  different 
phases  indicating  a  potential  source  of  bias  in  the  detennination  of  m(,  station  correc¬ 
tions  using  different  phases.  The  bias  is  expected  since  each  station  does  not  always 
detect  each  of  the  "a",  "b",  and  "max"  phase  for  an  event.  There  is  less  scattering  in 
the  station  effects  as  observed  between  the  "a"  and  "b"  phase  (Figure  6a)  compared  to 
those  between  the  "a"  and  "max"  phase  (Figure  6b),  indicating  that  there  is  a  greater 
common  number  for  "a"  and  ”b”  phase  observations  for  each  event  by  each  station 
than  of  "a"  and  "max"  phase  observations.  Figure  7  plots  the  LS-GLM  versus  the 
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Table  4.  iVILE-GLM88  Station  Corrections 


Station 

MLE-GLM 

correction 

a 

LS-GLM 

correction 

a 

#  of  observations 
signal  noise  clipped 

aae 

-0.190 

0.027 

-0.171 

0.033 

66 

71 

16 

aam 

0.279 

0.024 

0.239 

0.024 

125 

52 

6 

ade 

-0.005 

0.051 

0.128 

0.066 

16 

25 

0 

afi 

-0.148 

0.034 

0.007 

0.051 

27 

65 

0 

aku 

-0.135 

0.029 

-0.038 

0.034 

62 

63 

0 

alq 

0.030 

0.024 

-0.099 

0.022 

140 

22 

19 

alt 

-0.160 

0.190 

-0.209 

0.153 

3 

0 

0 

anp 

-0.376 

0.040 

-0.118 

0.066 

16 

53 

0 

ant 

0.056 

0.034 

0.097 

0.041 

41 

49 

2 

aqu 

-0.093 

0.030 

-0.202 

0.033 

66 

44 

14 

are 

0.104 

0.030 

0.038 

0.029 

83 

40 

0 

ath 

0.356 

0.134 

0.253 

0.119 

5 

1 

0 

atl 

0.150 

0.030 

0.067 

0.028 

90 

26 

2 

atu 

0.159 

0.025 

0.147 

0.027 

94 

66 

14 

bag 

-0.020 

0.025 

0.026 

0.025 

113 

52 

8 

bdf 

-0.011 

0.095 

-0.060 

0.084 

10 

2 

0 

bee 

-0.126 

0.029 

-0.029 

0.041 

41 

85 

3 

bhp 

-0.181 

0.032 

0.038 

0.049 

30 

75 

0 

bks 

0.083 

0.025 

0.080 

0.024 

124 

50 

1 

bla 

0.155 

0.024 

0.047 

0.023 

137 

39 

12 

bog 

0.222 

0.027 

0.336 

0.032 

69 

78 

0 

boz 

0.259 

0.051 

-0.009 

0.046 

34 

3 

4 

bul 

0.094 

0.026 

0.021 

0.024 

121 

28 

8 

car 

0.153 

0.026 

0.127 

0.028 

92 

59 

7 

ccg 

-0.125 

0.329 

-0.657 

0.266 

1 

0 

0 

chg 

-0.210 

0.031 

-0.368 

0.032 

71 

16 

27 

cmc 

-0.213 

0.040 

-0.287 

0.041 

42 

26 

0 

col 

0.069 

0.019 

-0.040 

0.018 

229 

45 

25 

cop 

0.161 

0.025 

0.193 

0.034 

61 

94 

14 

cor 

0.089 

0.030 

0.090 

0.032 

71 

47 

3 

eta 

0.235 

0.042 

0.123 

0.039 

46 

12 

4 

dag 

-0.075 

0.052 

-0.132 

0.054 

24 

15 

1 

dal 

0.191 

0.049 

0.158 

0.064 

17 

24 

4 

dav 

-0.246 

0.037 

0.032 

0.055 

23 

58 

0 

dug 

-0.002 

0.025 

-0.148 

0.022 

144 

16 

18 

eil 

0.062 

0.048 

-0.229 

0.063 

18 

3 

27 

ept 

0.003 

0.057 

-0.117 

0.049 

29 

2 

2 

esk 

0.104 

0.027 

0.209 

0.031 

75 

75 

2 

flo 

0.080 

0.033 

-0.133 

0.031 

72 

16 

9 

fvm 

-0.010 

0.046 

-0.110 

0.040 

44 

8 

0 
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Table  4.  (cont’d) 


Station 

MLE-GLM 

correction 

a 

LS-GLM 

correction 

a 

#  of  observations 
signal  noise  clipped 

gdh 

-0.167 

0.021 

-0.151 

0.023 

138 

110 

1 

geo 

0.016 

0.027 

0.065 

0.030 

81 

64 

2 

gie 

-0.195 

0.044 

0.041 

0.089 

9 

46 

0 

gol 

-0.264 

0.026 

-0.340 

0.023 

130 

20 

9 

grm 

0.028 

0.085 

1.029 

0.266 

1 

14 

0 

gsc 

0.065 

0.032 

-0.039 

0.031 

73 

17 

16 

gua 

-0.247 

0.022 

-0.128 

0.030 

77 

149 

0 

hkc 

-0.107 

0.028 

-0.040 

0.030 

79 

55 

0 

hlw 

-0.009 

0.033 

-0.254 

0.041 

43 

25 

29 

hnr 

0.228 

0.033 

0.380 

0.048 

31 

67 

0 

hon 

0.043 

0.085 

0.132 

0.108 

6 

9 

0 

ist 

0.189 

0.023 

0.162 

0.028 

88 

94 

25 

jet 

0.157 

0.035 

0.017 

0.035 

59 

4 

24 

jer 

0.057 

0.029 

0.046 

0.032 

68 

39 

20 

kbl 

0.146 

0.043 

-0.208 

0.071 

14 

0 

46 

kbs 

-0.154 

0.039 

-0.152 

0.041 

42 

30 

0 

kev 

-0.141 

0.024 

-0.120 

0.026 

104 

87 

4 

kip 

0.118 

0.023 

0.255 

0.029 

84 

118 

0 

kod 

0.040 

0.028 

-0.023 

0.030 

81 

29 

25 

kon 

0.089 

0.022 

-0.011 

0.026 

108 

61 

63 

krk 

-0.348 

0.095 

-0.330 

0.119 

5 

7 

0 

ktg 

-0.273 

0.029 

-0.199 

0.035 

58 

66 

1 

lah 

0.419 

0.066 

0.555 

0.119 

5 

17 

3 

lem 

-0.507 

0.031 

-0.338 

0.037 

51 

59 

0 

Ion 

-0.047 

0.024 

-0.138 

0.023 

133 

38 

21 

lor 

0.219 

0.035 

0.069 

0.033 

63 

8 

16 

Ipa 

0.419 

0.033 

0.868 

0.094 

8 

91 

0 

Ipb 

-0.042 

0.033 

-0.116 

0.035 

58 

39 

3 

Ips 

-0.069 

0.037 

-0.108 

0.038 

50 

27 

3 

lub 

0.210 

0.039 

0.174 

0.042 

40 

30 

3 

mal 

0.052 

0.030 

0.005 

0.030 

78 

37 

5 

man 

0.343 

0.060 

0.288 

0.058 

21 

8 

1 

mat 

-0.083 

0.024 

-0.136 

0.023 

130 

40 

18 

mds 

-0.087 

0.046 

-0.179 

0.046 

33 

19 

0 

mhi 

0.356 

0.110 

0.291 

0.119 

5 

2 

2 

mnn 

0.087 

0.076 

-0.170 

0.084 

10 

7 

2 

msh 

0.170 

0.042 

0.221 

0.046 

34 

19 

9 

mso 

0.059 

0.044 

-0.071 

0.039 

46 

7 

2 

mun 

-0.038 

0.043 

0.087 

0  050 

28 

31 

0 

nai 

-0.065 

0.028 

-0.110 

0.028 

93 

36 

9 

nat 

0.067 

0.045 

0.138 

0.051 

27 

27 

0 
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Table  4.  (cont’d) 


station 

MLE-GLM 

correction 

a 

LS-GLM 

correction 

CT 

#  of  observations 
signal  noise  clipped 

ndi 

0.108 

0.028 

0.028 

0.027 

100 

23 

20 

nha 

-0.214 

0.134 

-0.168 

0.133 

4 

2 

0 

nil 

-0.010 

0.053 

-0.076 

0.071 

14 

6 

18 

nna 

-0.163 

0.032 

-0.093 

0.039 

47 

57 

0 

nor 

-0.304 

0.032 

-0.247 

0.034 

62 

45 

0 

nur 

-0.017 

0.026 

-0.013 

0.030 

77 

77 

3 

ogd 

-0.090 

0.025 

-0.125 

0.024 

124 

42 

6 

oxf 

0.353 

0.030 

0.205 

0.028 

91 

16 

17 

pda 

0.074 

0.029 

0.195 

0.048 

31 

91 

3 

pel 

0.052 

0.039 

0.077 

0.049 

29 

39 

3 

pmg 

-0.082 

0.032 

-0.057 

0.034 

62 

45 

0 

pKX) 

0.006 

0.026 

-0.070 

0.025 

no 

33 

23 

pre 

-0.020 

0.038 

-0.118 

0.038 

49 

28 

0 

pto 

-0.127 

0.029 

-0.147 

0.030 

76 

45 

5 

que 

-0.410 

0.031 

-0.602 

0.036 

56 

22 

33 

qui 

0.002 

0.038 

0.188 

0.089 

9 

67 

0 

rab 

-0.152 

0.028 

0.072 

0.043 

38 

104 

0 

rar 

-0.106 

0.050 

0.080 

0.077 

12 

32 

0 

red 

0.484 

0.050 

0.373 

0.051 

27 

14 

3 

riv 

0.348 

0.059 

0.595 

0.089 

9 

22 

0 

sba 

-0.632 

0.088 

-0.320 

0.188 

2 

12 

0 

sep 

0.056 

0.022 

-0.009 

0.022 

148 

55 

18 

sdb 

0.070 

0.037 

-0.005 

0.036 

56 

15 

9 

seo 

-0.016 

0.026 

0.054 

0.029 

82 

65 

11 

sha 

0.343 

0.028 

0.358 

0.030 

76 

65 

0 

shi 

0.308 

0.032 

0.127 

0.032 

67 

11 

25 

shk 

-0.293 

0.032 

-0.069 

0.042 

40 

68 

0 

shl 

0.067 

0.032 

-0.137 

0.036 

55 

15 

34 

sjg 

-0.251 

0.024 

-0.244 

0.023 

129 

57 

0 

sna 

0.105 

0.080 

0.067 

0.108 

6 

11 

0 

sng 

-0.083 

0.043 

-0.035 

0.045 

35 

24 

1 

spa 

-0.759 

0.074 

-0.761 

0.074 

13 

7 

0 

stu 

0.055 

0.021 

0.037 

0.022 

142 

89 

20 

tab 

0.172 

0.034 

0.199 

0.036 

56 

37 

3 

tau 

-0.119 

0.065 

-0.067 

0.077 

12 

14 

0 

tol 

0.235 

0.026 

0.189 

0.028 

91 

42 

22 

tri 

-0.159 

0.023 

-0.209 

0.026 

105 

81 

17 

tm 

0.098 

0.024 

0.103 

0.025 

112 

70 

1 

tuc 

0.075 

0.037 

-0.062 

0.035 

57 

3 

21 

ume 

0.157 

0.027 

0.101 

0.026 

102 

48 

1 

unm 

-0.244 

0.067 

-0.026 

0.084 

10 

13 

1 
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Table  4.  (cont’d) 


Station 

MLE-GLM 

correction 

a 

LS-GLM 

correction 

a 

#  of  observations 
signal  noise  clipped 

upa 

-0.260 

0.190 

-0.406 

0.188 

2 

1 

0 

val 

0.009 

0.022 

0.053 

0.026 

104 

108 

12 

wel 

0.131 

0.074 

0.393 

0.094 

8 

12 

0 

wes 

-0.106 

0.022 

-0.043 

0.024 

124 

97 

6 

win 

-0.193 

0.049 

-0.149 

0.057 

22 

24 

0 
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STATION  EFFECTS  (censoring) 


Figure  6a.  MLE-GLM  station  effects  derived  using  the  "a"  phase  plotted  versus  those 
for  the  "b"  phase.  The  bias  between  the  two  estimates  is  up  to  0.2  magnitude  units. 
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Figure  6b.  MLE-GLM  station  effects  derived  using  the  "a"  phase  plotted  versus  those 
for  the  "max"  phase.  The  bias  between  the  two  estimations  is  much  higher,  being  up  to 
0.4  magnitude  units. 
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STAIION  EFFECTS  (censoring) 


Figure  6c.  MLE-GLM  station  effects  derived  using  the  "b"  phase  plotted  versus  those 
for  the  "max"  phase.  The  bias  between  the  two  estimates  is  up  to  0.2  magnitude  units. 
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STATION  EFFECTS 


Figure  7.  LS-GLM  station  effects  versus  MLE-GLM  station  effects  for  127  WWSSN 
stations.  Note  that  for  most  stations  the  difference  is  skewed  and  therefore  potential 
for  bias  exists  when  the  censoring  information  is  ignored. 
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MLE-GLM  station  effects  derived  for  all  magnitude  models.  The  net  difference 
between  the  two  data  sets  is  by  definition  zero,  but  the  scatter  in  the  results  does  not 
allow  one  to  infer  such  from  the  plot.  The  possible  departure  of  the  slope  from  unity 
in  fitting  the  data  in  Figure  7  indicates  the  degree  of  bias  in  the  LS-GLM  station 
corrections  when  comparing  to  the  MLE-GLM  model  which  requires  no  apriori  sta¬ 
tion  corrections  prior  to  inversion.  Based  on  Ringdal’s  (1986)  LS-GLM  station  terms, 
Veith  et  al.  (1987)  computed  detection  thresholds  that  are  undistinguishable  at  noisy 
and  quiet  stations.  The  bias  in  LS-GLM  has  been  demonstrated  by  Lilwall  (1986)  to 
be  attributed  to  the  ambiguities  in  keeping  the  detection  threshold  constant. 
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3.5  Error  Analysis 

The  uncertainty  associated  with  each  event  and  station  is  computed  with  two 
methods.  The  standard  deviation  a  of  the  perturbing  noise  returned  by  the  EM  algo¬ 
rithm  is  believed  to  be  a  good  measure  of  the  uncertainty  in  any  single  observation. 
Jih  et  al.  (1988)  propose  to  scale  this  a  with  the  number  of  associated  observations  to 
infer  the  precision  of  each  estimated  parameter.  The  second  method  is  by  utilizing  the 
bootstrap  method  (Efron,  1979,  1981),  which  has  been  extensively  used  in  some  recent 
network  m^  estimation  studies  (McLaughlin  etal.,  1986a,  1986b,  1988).  Jih  et  al. 
(1988)  propose  to  enlarge  the  pool  of  regular  residuals  with  the  generalized  (i.e.  cen¬ 
sored)  residuals  for  resampling  to  account  for  the  fact  that  all  the  paths  containing 
clipped  or  noisy  data  to  remain  so  throughout  the  resampling  iteration. 

The  m^  error  estimates  are  plotted  in  Figure  8  using  the  bootstrap  and  EM  for 
censored  data  only.  Figure  9  shows  the  m^  error  estimates  obtained  using  censoring 
and  no-censoring  data,  respectively.  The  EM  nij,  error  estimates  using  censoring  data, 
as  shown  in  Figure  9a,  are  generally  larger  than  the  no-censoring  data.  The  same  is 
true  for  the  estimates  using  bootstrap  as  shown  in  Figure  9b.  The  bootstrap  method, 
requiring  intensive  computing  time,  does  not  yield  correct  mb  estimates  and  is  used 
only  for  error  estimation  (e.g.  Blandford  et  al.  1984,  McLaughlin,  1988).  Taking  into 
account  the  fact  that  EM  estimation  indeed  yields  the  correct  mb  values  together  with 
its  standard  errors,  we  prefer  to  use  this  method  over  bootstrapping.  The  resampling 
scheme  used  in  the  bootstrap  estimation  has  the  drawback  of  assuming  that  all  events, 
regardless  of  its  size,  may  be  censored  with  equal  chance.  This  may  only  be  true  when 
events  of  essentially  the  same  mb  level  are  used  in  the  MLE-GLM.  This  assumption 
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Figure  9a.  LS-GLM  mj,  error  estimates  versus  MLE-GLM  mj,  error  estimates  obtained 
from  EM.  The  LS-GLM  mj,  errors  computed  with  no-censoring  information  are  higher 
than  those  with  censoring  information. 
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Figure  9b.  LS-GLM  nib  error  estimates  versus  MLE-GLM  mb  error  estimates  obtained 
using  bootstrap.  The  LS-GLM  mb  error  estimates  computed  with  no-censoring  infor¬ 
mation  are  higher  than  those  with  censoring  information. 
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may  give  rise  to  some  of  the  discrepancies  in  the  mb  error  estimation  using  the  EM 
and  bootstrap  approaches.  The  MLE-GLM  mb  error  estimates  for  censoring  and  no- 
censoring  data  are  plotted  in  Figure  10  for  all  three  phases  separately.  The  larger 
degree  of  deviation  from  a  unit  slope  for  the  error  of  mb(a)  than  for  the  other  two 
phases  may  be  attributed  to  the  smaller  number  of  observations  for  phase  "a". 

The  error  of  the  station  effects  is  also  computed.  As  in  the  case  of  mb  error  esti¬ 
mates,  the  MLE-GLM  error  estimates  of  the  station  effects  using  censoring  data  are 
smaller  than  those  for  GLM  using  no-censoring  information,  as  shown  in  Figure  11a. 
Using  the  station  error  estimates  from  bootstrap,  the  deviation  from  those  obtained 
using  no-censoring  information  becomes  less  than  the  MLE-GLM  case,  as  shown  in 
Figure  lib. 
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m^(a)  error 


Figure  10a.  LS-GLM  ntj,  error  estimates  versus  MLE-GLM  mj,  error  estimates  obtained 
from  EM  using  just  the  "a"  phase  data.  The  LS-GLM  m^  error  estimates  computed 
with  no-censoring  information  are  higher  than  those  with  censoring  information. 
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Figure  10b.  LS-GLM  nib  error  estimates  versus  MLE-GLM  mb  error  estimates 
obtained  from  EM  using  just  the  "b"  phase  data.  The  LS-GLM  mb  error  estimates 
computed  with  no-censoring  information  are  higher  than  those  with  censoring  informa¬ 
tion  but  show  smaller  bias  than  those  for  phase  "a". 


Tcledyne  Gcotech 


57 


June  1988 


censoring  (EM) 


Final  Report  Magnitude  Yield  Estimation  TGAL-87-05 


m^(max)  error 

0.20 


0.15 


0.10 


0.05 


0.00 

0.00  0.10  0.20 

no  — censoring 

Figure  10c.  LS-GLM  mb  error  estimates  versus  MLE-GLM  mb  error  estimates 
obtained  from  EM  using  just  the  "max"  phase  data.  The  LS-GLM  mb  error  estimates 
computed  with  no-censoring  information  are  higher  than  those  with  censoring  informa¬ 
tion  but  show  smaller  bias  than  those  for  phase  "a". 
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Figure  11a.  LS-GLM  site  term  error  estimates  versus  MLE-GLM  site  term  error  esti¬ 
mates  using  EM.  The  error  estimates  for  the  no-censoring  case  are  larger  than  those 
for  the  censoring  case  in  agreement  with  the  m(,  error  observations. 
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Figure  11b.  MLE-GLM  site  term  error  estimates  obtained  using  bootstrap  are  plotted 
against  those  obtained  using  EM.  The  error  estimates  using  the  two  different  schemes 
are  in  quite  close  agreement  to  each  other. 
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4.0  PHASE  DIFFERENTIAL  nib's 

P  and  P  coda  arrivals  arc  most  likely  affected  by  anomalous  arrivals  due  to 
scattering,  surface  multiples,  and  cratering.  The  first  motion  P-wave  amplitude  meas¬ 
ured  as  the  "a"  phase  is  perhaps  the  least  affected  by  these  later  arrivals  compared  to 
the  "b"  or  "max"  phase.  Therefore,  comparing  the  "a"  phase  to  the  "b"  or  "max"  phase 
may  provide  a  calibration  on  the  near-source  effects  of  a  nuclear  explosion.  By  using 
the  MLE-GLM  mb  estimates  from  the  phase-ratios  of  these  phases,  characteristics  of 
the  lateral  variability  of  the  near-source  structure  will  be  examined. 

The  log  amplitude  ratios  between  the  "a"  and  "max"  phase  have  been  studied  to 
discriminate  between  contained  and  cratering  events  at  Shagan  River  and  Degelen  by 
McLaughlin  et  al.  (1985,  1986a).  Their  findings  were  that  the  log  amplitude  ratios  of 
the  "max"  phase  over  "a"  phase  for  cratering  explosions  are  smaller  than  the  contained 
ones.  The  coda  energy  level  of  the  teleseismic  P-wave  decreases  dramatically  for 
cratering  shots.  An  example  of  such  is  the  Jan  15,  1965  presumed  cratering  explosion 
at  Shagan  River  (Rodean,  1979).  Using  this  event  as  a  calibration,  McLaughlin  et  al. 
(1986a)  have  found  two  other  Shagan  River  events  which  have  similar  amplitude  ratios 
approaching  that  of  the  Jan  15,  1965  event. 

In  this  study,  the  MLE-GLM  mb  estimates  for  NTS  and  Novaya  Zemlya  explo¬ 
sions  are  analyzed  along  with  the  E.  Kazakhstan  explosions.  The  MLE-GLM  mb’s 
estimates  for  the  different  test  sites  are  used  to  compute  the  mb(max)-mb(a)  and 
mb(max)-mb(b)  for  all  events  of  interest  and  are  analyzed  in  the  context  of  each  test 
site.  A  self-consistent  set  of  station  terms  derived  in  the  MLE-GLM  is  used  to  correct 
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for  the  site  effects  in  the  maximum-likelihood  estimation. 

The  comparisons  between  the  different  phases  made  in  Figure  3  for  all  events 
show  that  the  m^  estimates  for  the  "a"  phase  are  on  the  average  0.5  magnitude  unit 
less  than  those  for  the  "max"  phase.  The  scatter  in  the  data  indicates  that  biases  exist 
among  the  different  events,  which  warrants  a  more  in-depth  analysis.  The  Sm^’s 
between  "a"  and  "max"  and  "b"  and  "max"  phases  are  plotted  against  mb(max)  in  Fig¬ 
ures  12a  and  12b  for  Shagan  and  Degelen,  respectively.  Considering  only  those  events 
recorded  by  more  than  10  stations  for  the  "a"  phase,  the  5mi,(max-a)’s  range  between 
0.4  to  0.68  for  Shagan  and  0.45  to  0.54  for  Degelen  events.  These  values  are  within 
the  range  presented  in  McLaughlin  et  al.  (1985).  The  smaller  range  in  Smb’s  obtained 
for  the  Degelen  events  compared  to  Shagan  may  imply  that  the  structure  beneath 
Degelen  is  more  nearly  homogeneous,  but  the  sample  is  too  small  to  be  conclusive. 

The  Jan  15,  1965  (event  29  in  Table  2  and  Figure  12a)  presumed  cratering  event 
shows  a  lower  than  average  6mb(max-a)  of  0.39,  in  agreement  with  the  observation  by 
McLaughlin  et  a/.  (1985).  But  further  analysis  of  this  event  indicates  that  its  5mb(max- 
b)  does  not  deviate  statistically  from  the  rest  of  the  population,  possibly  indicating  that 
the  "a"  phase  is  less  affected  by  the  surface  reflections  than  the  "b"  and  "max"  phases. 
The  Sep  15,  1978  (event  38)  explosion  which  shows  an  anomalous  low  log(max/a)  in 
McLaughlin  et  al.  (1985)  does  not  appear  to  be  such  in  our  analysis.  This  event  has  a 
6mb(max-a)  of  0.47  which  is  not  statistically  different  from  the  rest  of  the  5mb(max- 
a)’s  taking  into  account  2a  error  of  about  0.1 1  from  Table  3.  But  in  fact,  the  5mb  for 
this  event  is  even  higher  than  two  other  events,  Jul  23,  1973  (event  32)  and  Jun  23, 
1979  (event  39),  that  were  not  studied  in  McLaughlin  et  al.  ’s  (1985).  The  scatter  in 
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Figure  12a.  m5(max)’s  plotted  against  6mb’s  for  Shagan.  The  mb’s  are  maximum 
likelihood  estimates  with  the  censoring  information.  The  solid  circles  are  mb(max)- 
mb(a)  and  the  open  circles  are  mb(max)-mb(b)  for  the  corresponding  event.  The  events 
are  listed  in  Table  2. 
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Figure  12b.  ni(,(max)’s  plotted  against  Sm^’s  for  Degelen.  The  nib’s  are  maximum 
likelihood  estimates  with  the  censoring  information.  The  solid  circles  are  mb(max)- 
m(,(a)  and  the  open  circles  are  mb(max)-m(,(b)  for  the  corresponding  event.  The  events 
are  listed  in  Table  2.  The  hatched  circles  are  mb(max)-mb(a)  for  which  there  are  less 
than  10  "a"  phase  signal  readings. 
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the  data  nonetheless  prevents  us  to  argue  with  high  confidence  that  the  low  6mb(max- 
a)  of  event  29  is  directly  related  to  the  cratering  effect.  Moreover,  the  only  physical 
interpretation  of  such  a  relationship  is  possibly  that  of  an  absence  of  surface  reflected 
multiples  causing  a  reduction  in  the  "max"  and  "b"  amplitudes.  Such  an  effect  was 
however  not  observed  in  a  comparison  of  regional  from  cratering  and  tamjjed 
explosions  (Gupta  et  al.,  1988). 

For  the  NTS  events,  most  of  the  explosions  with  mb(max)  of  less  than  5.0  are  not 
well  detected  by  the  WWSSN  network  and  are  not  used  in  the  analysis.  The  Sm^’s 
are  plotted  in  Figure  12c.  All  events  with  m|j(max)  of  less  than  5.5  have  low 
5mb(max-a)’s  between  0.3  to  0.43.  The  rest  of  the  larger  events  have  5mb(max-a)’s  of 
over  0.5,  except  for  Handley  (event  11).  The  anomalous  low  magnitude  ratios  for  the 
smaller  events  (mb(max)  between  5.0  to  5.5)  may  be  attributed  to  low  detection  levels, 
and  they  are  therefore  not  emphasized  in  the  analysis.  Handley  (mb(max)  6.48),  which 
is  well  recorded  by  the  network  with  41  signal  readings,  exhibits  an  anomalously  large 
"a"  amplitude  which  indicates  that  there  may  be  strong  heterogeneity  in  its  source 
function  and  source  medium  beneath  this  test  site  compared  to  the  rest  of  the  NTS  test 
sites. 

The  Novaya  Zemlya  events,  divided  into  north  and  south,  arc  plotted  in  Figure 
12d.  The  presumed  double  event  of  Oct  18,  1975  (event  55)  (Hurley,  1977;  Chan 
et  al.  1988)  stands  out  as  having  large  coda  multiples,  which  contribute  to  the 
anomalous  large  "max"  phase  compared  to  the  "a"  phase.  Most  of  the  north  Novaya 
Zemlya  events  have  5mb’s  of  between  0.35  to  0.48,  which  are  much  less  than  the  other 
test  sites,  indicating  that  the  source  structure  there  may  be  quite  different  from  the 
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Figure  12c.  mb(max)’s  plotUid  against  5mb’s  for  NTS.  The  mb’s  are  maximum  likeli¬ 
hood  estimates  with  the  censoring  information.  The  solid  circles  are  mb(max)-mb(a) 
and  the  open  circles  are  mb(max)-mb(b)  for  the  corresponding  event.  The  events  are 
listed  in  Table  2.  The  hatched  circles  are  mb(max)-mb(a)  for  which  there  are  less  than 
10  "a"  phase  signal  readings. 
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Figure  12d.  mb(inax)’s  plotted  against  Sm^’s  for  Novaya  Zemlya.  The  nib’s  are  max¬ 
imum  likelihood  estimates  with  the  censoring  information.  The  solid  and  open  circles 
are  mb(max)-mb(a)  and  mb(max)-mb(b),  respectively,  for  corresponding  events  in  the 
northern  site.  The  solid  and  open  squares  are  m|,(max)-mi,(a)  and  mb(max)-mb(b), 
respectively,  for  corresponding  events  in  the  southern  site.  The  hatched  circles  and 
squares  are  mb(max)-mb(a)  for  which  there  are  less  than  10  "a"  phase  signal  readings. 
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Others.  There  is  a  slight  indication  that  the  structure  beneath  the  southern  and  northern 
sites  may  be  different  judging  from  their  5m|,’s  levels,  but  more  data  from  this  test  site 


need  to  be  analyzed  to  comment  further  on  this  issue. 
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5.0  SINGLE  TEST  SITE  GLM’s  VERSUS  WORLD-WIDE  GLM’s 

As  was  stated  earlier,  Longshot  was  not  an  adequate  calibration  event  for  Milrow 
and  Cannikin  because  the  composition  of  WWSSN  stations  changed  for  the  three 
events  and  because  of  the  clipping  for  Cannikin  and  Milrow.  In  order  to  demonstrate 
this,  we  list  in  Table  5  the  magnitudes  for  these  three  shots  as  determined  from  an 
LS-GLM  and  MLE-GLM  using  only  these  three  events,  compared  to  the  results  of  the 
full  MLE-GLM  for  m^Cmax).  The  results  are  striking  in  that  Cannikin  decreases  by  up 
to  0.105  magnitude  units  from  a  single  test-site  determination  to  the  full  determination, 
whereas  Longshot  increases  by  0.012. 


Table  5.  Comparison  of  Single  Test  Site  Only  GLM  with  MLE-GLM88 

Amchitka  Test  Site 


Amchitka  Only  World-Wide 


Event 

LS-GLM 

MLE-GLM 

LS-GLM 

MLE-GLM 

Smle-cjlm 

Cannikin 

6.984 

7.021 

6.832 

6.916 

0.105 

Milrow 

6.565 

6.583 

6.433 

6.497 

0.086 

Longshot 

5.821 

5.811 

5.854 

5.823 

-0.012 

This  comparison  points  out  that  the  common  practice  of  choosing  only  a  few  cali¬ 
bration  events  at  a  single  test  site  can  and  does  lead  to  systematic  errors.  How  does 
this  come  about?  When  we  examine  the  Amchitka  data  set,  we  see  that  it  contains  84 
stations  and  that  the  RMS  station  effect  is  0.37  magnitude  units.  If  we  were  to  assume 
that  station  effects  are  normal  random  variables  with  a  standard  deviation  of  0.37  and 
zero  mean,  then  we  would  expect  that  the  mean  of  84  stations  drawn  from  the  popula¬ 
tion  would  have  a  5%  probability  of  being  either  greater  than  0.08  or  less  than  -0.08. 
In  actuality  there  are  only  32  common  WWSSN  stations  (20°  <  A  <  QS”!  between 
Longshot  and  the  two  larger  shots  Milrow  and  Cannikin.  The  95%  confidence  bound 
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on  the  bias  from  these  32  stations  is  ±  0.13  magnitude  units. 

It  is  important  to  remember  that  the  constraint  equation  that  the  station  effects 
N 

should  sum  to  zero,  =  0,  has  an  uncertainty  associated  with  it.  In  the  case  of 

k=l 

the  Amchitka  explosions,  this  potential  bias  appears  to  be  more  acute  than  at  other  test 
sites  since  the  network  configuration  had  changed  substantially  between  these  explo¬ 
sions  in  the  late  1960’s. 

Of  course,  in  reality  some  of  these  biases  are  caused  by  deterministic  effects  that 
one  can  anticipate.  The  leading  hypothesis  for  this  Amchitka  "bias"  is  the  focusing- 
defocusing  of  the  P-wave  by  the  descending  slab  beneath  the  Aleutian  island  arc.  Such 
effects  could  be  predicted  deterministically  and  corrections  could  be  applied,  as  was 
done  by  Cormier  (1986)  using  ray-tracing  methods  applied  to  a  3-D  model  of  the 
source  region  obtained  from  block  inversion  of  traveltime  residuals. 

We  have  performed  the  same  analysis  for  9  NTS  events,  16  Novaya  Zemlya 
events,  and  11  Tuomotu  events  with  the  results  that  the  event  magnitudes  had  RMS 
changes  of  0.038,  0.072,  and  0.025  for  NTS,  Novaya  Zemlya,  and  Tuomotu  respec¬ 
tively  (see  Table  6).  None  of  these  test  sites  demonstrated  substantial  changes 
approaching  0.1  magnitude  units.  In  each  case,  the  data  were  analyzed  using  an 
MLE-GLM  for  the  test  site  alone,  and  then  the  individual  event  magnitudes  were  com¬ 
pared  to  an  MLE-GLM  with  a  full  complement  of  106  explosions.  However,  when  6 
southern  Novaya  events  were  used  in  an  MLE-GLM,  the  average  event  decreased  0.05 
magnitude  units  with  an  rms  shift  of  0.077  magnitude  units.  This  is  significant  in  light 
of  the  formal  standard  errors  of  the  individual  event  magnitudes  of  about  0.02  units. 
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This  implies  that  if  magnitude  calibrations  are  to  be  performed  for  a  single  test  site, 
then  uncertainties  approaching  0.1  magnitude  units  should  be  attributed  to  the  final 
results  for  purposes  of  inter-test  site  comparison. 

It  is  interesting  to  note  that  the  RMS  station  effect  and  the  residual  RMS  error  are 
roughly  in  the  same  proportion  for  each  test  site  that  was  examined  separately.  How¬ 
ever,  there  may  be  significant  differences  between  some  test  sites  in  the  total  variance. 
Table  6  lists  the  number  of  events,  the  number  of  stations,  the  degrees  of  freedom,  the 
RMS  station  effect,  the  RMS  residual  error,  and  Oiotaj  =  V (RMS-STA)^-t-(RMS-RES)^ 
for  each  major  test-site  specific  MLE-GLM  and  an  MLE-GLM  based  on  111  explo¬ 
sions.  We  see  that  the  station  effects  are  between  73%  and  90%  of  the  original  magni¬ 
tude  variance.  The  RMS  residual  is  smaller  for  smaller  sets  of  explosions  (fewer 
degrees  of  freedom)  and  for  sets  of  events  located  in  smaller  test  sites.  For  com¬ 
parison,  an  MLE-GLM  for  111  explosions  and  127  stations  had  2598  degrees  of  free¬ 
dom,  an  RMS  station  effect  of  0.219,  and  an  RMS  residual  of  0.247.  For  the  world¬ 
wide  MLE-GLM  the  station  effects  were  on  average  44%  of  the  total  original  variance. 
As  events  are  added  over  a  greater  area  of  the  globe,  the  station  effects  are  a  smaller 
proportion  of  the  total  variance,  tending  to  about  40%.  For  example,  the  station 
effects  for  northern  and  southern  Novaya  Zemlya  test  sites  treated  separately  are  only 
53%  correlated;  yet  the  two  testing  areas  are  separated  by  only  800  km.  Amchitka 
(AMC)  and  Northern  Novaya  Zemlya  (NNZ)  have  the  largest  total  variance,  and  NTS 
has  the  smallest. 
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Table  6.  Comparison  of  Major  Test  Site  Specific  mb(max)  MLE-GLM’s 

Test  Site 

#  Events 

#  Stations 

#  Deg.  Freed. 

RSM  STA 

RMS  RES. 

®total 

NTS 

9 

64 

231 

0.246 

0.150 

0.288 

SNZ 

6 

43 

45 

0.273 

0.138 

0.306 

TUO 

11 

53 

223 

0.286 

0.122 

0.311 

NNZ 

10 

75 

217 

0.357 

0.140 

0.383 

NZ 

16 

76 

303 

0.348 

0.169 

0.387 

AMC 

3 

84 

82 

0.370 

0.124 

0.390 

1 1  lEXP 

111 

127 

2598 

0.219 

0.247 

0.339 
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6.0  INTER-EVENT  CORRELATION 

The  procedure  by  Blandford  et  al.  (1984)  for  perturbing  the  MLE-GLM  solution 
with  an  a  priori  correlation  matrix  is  used  by  McLaughlin  et  al.  (1986a)  to  construct 
GLM86  for  a  world-wide  set  of  53  events  included  in  the  MLE-GLM88  set.  They 
used  an  event-event  correlation  model,  so  that  the  station  corrections  would  not  be 
overly  influenced  by  a  large  number  of  events  at  a  single  test  site.  The  model  assumes 
that  pairs  of  observations  at  a  specific  station  are  correlated  if  the  two  observations  are 
from  events  that  are  close  together.  The  correlation  function  is  parameterized  by  two 
exponential  decays.  If  two  events,  i  and  j,  are  separated  by  distance  Ajj,  then  the 
observations  at  any  station  for  these  two  events  are  correlated  by  the  amount 

Cij  =  Aq  (exp(-Aij/ai)  exp(-Aij/a2))/2. 

where  Aq  =  0.8,  aj  =  20.0  km,  and  -  ^00  km  in  this  study.  By  definition,  Cy  =  1  if 
the  events  are  the  same  event,  i  =  j. 

The  procedure  is  approximate  and  would  not  be  needed  if  the  events  were  evenly 
distributed  around  the  world  and  each  station  were  to  record  signals  from  many  events 
with  an  unbiased  geographic  distribution.  In  practice,  if  a  large  number  of  events  are 
selected  from  as  many  test  sites  as  possible,  the  results  of  the  perturbation  could  be 
minimal.  To  test  this,  the  mb(max)  magnitudes  for  111  explosions  were  analyzed 
using  an  MLE-GLM  with  and  without  the  event-event  correlation  correction. 

Plots  of  the  estimates  derived  with  and  without  the  event-event  correlation 
perturbation  are  shown  in  Figure  13  for  the  three  phases  separately.  From  Figure  13, 
the  only  events  that  change  by  more  than  0.075  units  are  the  8  Azgir  events  that  were 
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Figure  13a.  Event  magnitudes  for  phase  "a",  mj,(a),  from  MLE-GLM  solutions  assum¬ 
ing  event-event  correlation  versus  no  event-event  correlation. 
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m^,(b) 


Figure  13b.  Event  magnitudes  for  phase  "b",  m|,(b),  from  MLE-GLM  solutions  assum¬ 
ing  event-event  correlation  versus  no  event-event  correlation. 
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Figure  13c.  Event  magnitudes  for  phase  "max",  mj,(max),  from  MLE-GLM  solutions 
assuming  event-event  correlation  versus  no  event-event  correlation. 
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not  read  at  the  full  complement  of  stations.  These  events  were  read  at  a  small  number 
of  stations,  and  therefore  perturbations  to  the  small  set  of  station  corrections  can  intro¬ 
duce  statistically  significant  perturbations.  The  2  Azgir  events  that  were  read  at  the  full 
complement  of  stations  did  not  change  by  more  than  0.05.  Again  this  illustrates  the 
effects  of  the  statistics  of  offsets  in  magnitudes  from  networks  with  small  sample  size. 
The  m^  estimates  for  the  different  phases  obtained  with  event-event  correlation  are  also 
plotted  in  Figure  14.  A  similar  phase  bias  pattern  is  seen  here  as  for  those  in  the  EM 
case  (Figure  3).  The  station  effects  derived  with  and  without  event-event  correlation 
are  plotted  in  Figure  15.  The  scatter  of  the  station  effects  plotted  in  Figure  15a 
between  the  no-censoring  case  (LS-GLM)  and  the  no-censoring  event-event  correlation 
case  is  very  similar  to  that  between  the  EM  censoring  and  no-censoring  case  (Figure 
7).  The  station  effects  estimates  with  and  without  event-event  correlation  for  censor¬ 
ing  data  are  plotted  against  each  other  in  Figure  15b. 

The  resulting  event  magnitudes  differ  from  the  MLE-GLM  tables  presented  above 
(Table  2c)  without  the  correlation  correction  by  an  rms  amount  of  0.05  magnitude 
units.  The  station  effects  were  similarly  different  by  an  rms  difference  of  0.05  magni¬ 
tude  units.  The  differences  were  generally  small  and  comparable  to  the  differences 
between  the  MLE-GLM88  based  on  "a",  "b",  and  "max"  together.  Given  this  result, 
we  conclude  that  the  event-event  correlation  perturbation  is  not  required  for  MLE- 
GLM88  and  that  the  station  effects  are  not  significantly  improved  by  the  procedure. 
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Figure  14a.  m(,’s  estimates  assuming  event-event  correlation  for  the  "a"  phase  plotted 
against  the  "b"  phase.  The  bias  is  similar  to  that  observed  for  the  MLE-GLM  case. 


Teledyne  Geotech 


June  1988 


Final  Report 


Magnitude  Yield  Estimation 


TGAL-87.05 


m^(event  — correlation) 


4.0  5.0  6.0  7.0 

^,3(0) 


Figure  14b.  nib’s  estimates  assuming  event-event  correlation  for  the  "a"  phase  plotted 
against  the  "max"  phase.  The  bias  is  similar  to  that  observed  for  the  MLE-GLM  case. 
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Figure  15a.  Site  effects  for  MLE-GLM  solutions  assuming  event-event  correlation 
versus  LS-GLM  with  no-censoring  information  based  on  mj,’s  for  all  phases. 
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Figure  15b.  Site  effects  for  MLE-GLM  solutions  assuming  event-event  correlation 
versus  no  event-event  correlation  based  on  m^’s  for  all  phases. 
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7.0  CONCLUSIONS 

The  event  magnitudes  for  111  nuclear  explosions  at  over  eight  test  sites  world¬ 
wide  are  obtained  using  a  maximum-likelihood  estimation  and  constructed  in  a  general 
linear  model  (MLE-GLM88).  The  magnitude  estimates  and  station  corrections  are 
obtained  for  the  "a"  phase,  the  "b"  phase,  and  the  "max"  phase  of  the  short-period 
teleseismic  P-wave  using  data  from  the  WWSSN  network. 

•  Taking  into  account  the  station  censoring  information,  considerable  bias  is 
apparent  for  the  large  explosions  due  to  censoring  by  clipping,  and  for  small 
events  due  to  censoring  by  non-detection. 

•  When  censoring  information  is  considered,  the  m^  bias  is  0.28  for  mb(a)  vs 
mb(b),  0.19  for  m^Cb)  vs  m^Cmax),  and  0.48  for  m|j(a)  mb(max).  When  censor¬ 
ing  information  is  not  considered,  the  m^  biases  are  on  the  average  lower  from  a 
more  scattered  distribution. 

•  The  station  corrections  for  127  WWSSN  stations  will  be  useful  for  the  estimation 
of  seismic  magnitudes  for  other  events  not  included  in  MLE-GLM88.  It  would  be 
desirable  in  the  future  not  to  perurb  the  magnitudes  of  each  previously  deter¬ 
mined  event  as  new  events  are  analyzed.  In  this  regard,  the  WWSSN  network 
constitutes  a  basis  for  comparison. 

•  Error  estimations  were  performed  using  both  the  bootstrap  and  the  EM  algo¬ 
rithms.  The  mj,  errors  estimated  using  bootstrap  are  generally  larger  than  those 
obtained  from  EM.  The  error  estimates  taking  into  account  censoring  information 
are  generally  larger  than  those  obtained  without  considering  censoring 
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information. 

•  This  set  of  1 1 1  events  can  be  used  to  calibrate  smaller  networks  for  the  purposes 
of  test-site  specific  station  corrections.  Test  site  specific  station  corrections  yield 
more  precise  relative  magnitudes  for  a  test  site  but  are  susceptible  to  offsets 
between  test  sites.  However,  care  should  be  taken  in  deriving  station  corrections 
from  a  small  number  of  events  with  limited  geographic  extent. 

•  The  5mb’s  between  the  "max"  and  "a"  phases  for  a  presumed  cratering  event  may 
be  low  compared  to  the  contained  ones  due  to  the  weak  surface  reflections  associ¬ 
ated  with  a  cratered  explosion. 

•  The  6mb’s  between  the  "max"  and  "a"  phases  for  a  presumed  double  event  may 
be  larger  than  average  due  to  the  secondary  arrivals. 

•  The  potential  for  bias  in  determining  magnitudes  from  a  single  test-site  LS-GLM 
or  even  a  single  test  site  MLE-GLM  was  demonstrated  for  the  Amchitka  test  site. 
The  magnitudes  of  Cannikin  and  Milrow  may  be  biased  by  0. 1 1  magnitude  units 
high  due  to  the  geometry  of  the  WWSSN  network  if  a  reduction  is  performed  on 
the  Amchitka  test  site  only.  Station  corrections  for  a  world-wide  data  set  of 
explosions  reduce  the  magnitudes  of  the  largest  Amchitka  explosions  by  about 
0.12  magnitude  units. 

•  Because  the  "max"  phase  magnitudes  depend  so  strongly  on  the  non-linear  free- 
surface  interaction  (pP-nP),  the  "a"  and  "b"  phase  magnitudes  are  useful  in  estab¬ 
lishing  relative  sizes  of  events  with  less  concern  for  the  independent  estimation  of 
pP/P  relative  amplitudes  and  pP  delay  times  necessary  to  model  the  mb(max) 
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magnitude  for  a  given  test  site. 

•  Based  on  "a",  "b",  and  "max"  phase  magnitudes,  the  two  southern  Novaya  Zem- 
lya  events  of  02  November  1974  and  27  October  1973  are  0.1  and  0.2  magnitude 
units  (respectively)  larger  than  the  Amchitka  explosion  Cannikin.  All  three 
phases  are  consistent  in  this  regard,  and  the  conclusions  based  on  the  "a"  and  "b" 
phase  are  not  dependent  upon  the  interpretation  of  pP+P  interference  effects  on 
the  P  waveform.  The  two  northern  Novaya  Zemlya  events  of  12  September  1973 
and  14  October  1970  are  within  0.06  magnitude  units  of  the  same  size  as  Canni¬ 
kin  based  on  "a"  and  "b"  phase  magnitudes.  The  "max"  phase  magnitudes  of 
these  two  events  are  smaller  than  Cannikin.  Since  the  "max"  phase  magnitudes 
depend  on  P+pP  interference  effects,  the  "a"  and  "b"  phase  magnitudes  suggest 
that  these  two  events  are  the  same  size  as  Cannikin  and  that  the  lack  of  pP+P 
constructive  interference  is  responsible  for  the  smaller  "max"  phase  magnitudes  of 
the  largest  events  at  the  northern  Novaya  Zemlya  test  site. 

•  By  incorporating  inter-event  correlation  into  the  MLE-GLM88,  the  event  magni¬ 
tude  estimates  differ  by  an  rms  of  less  than  0.05  magnitude  units  and  is 
insignificant  compared  to  the  overall  error  of  the  model.  This  indicates  that  the 
inter-event  correlation  perturbation  may  not  be  necessary  when  a  large  database  is 
available. 
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XAMK 

mlglm  —  maximum  likelihood  general  linear  model 

SYNOPSIS 

migim  <  input 

DKSCRIPTION 

mIgIm  simultaneously  computes  maximum  likelihood  estimate  of  event  magnitudes  as  well  as  station 
effects  with  data  set  including  censored  information  and  inter-event  correlation.  In  fact,  it  is  applicable 
to  any  (general)  linear  models  of  the  form:  M(i,j)  »  E{i)  +  S(j)  +  n(i,j)  where  E()  and  S()  are  unk¬ 
nowns  to  be  estimated  and  M(,)  the  observed  data  matrix. 

It  uses  E-M  algorithms  of  Dempster  et  al.  (1977)  twice  to  solve  the  multiparameter  version  of  a  max¬ 
imum  likelihood  estimation  problem  originally  considered  by  Ringdal  in  1976.  In  the  first  loop,  the 
maximum  likelihood  estimates  are  computed  based  on  observations  without  censoring  information. 
This  result  is  merely  for  comparison  purpose  only.  While  in  the  second  loop,  program  recomputes  the 
estimates  with  whole  data  set  including  the  censoring  information  due  to  clipping  and  non-detection  as 
described  in  VSC-TR-82-12.  It  is  optional  to  adjust  the  result  with  inter-event  correlation  (TGAL-TR- 
83-06).  Efron’s  (1979,1981)  bootstrap  method  is  then  used  to  compute  the  uncertainty  in  the  estimates. 
Old  versions  of  this  code  have  been  utilized  extensively  in  the  past  four  years  (c.h  VSC-TR-82-12, 
TGAL-TR-83-06,  TGAL-86-01)  and  current  version  is  called  mlglm6  which  has  been  ported  to  and 
tested  on  various  types  of  computers  including  Celerity  12600,  VAX  11/780  and  SUN  3/160  at  Center 
for  Seismic  Studies  . 

A  simplified  version  of  the  code,  mlglm7,  is  also  available  in  which  we  always  bypass  the  incorporation 
of  event-event  correlation  function  with  the  linear  model.  This  is  due  to  the  urgent  need  of  processing 
of  large  database  (about  300  events)  which  is  far  beyond  the  capacity  of  migimfi. 

SAMPLE  INPUT  FILE 

Pmax,GLM_test,  9/1/1987 
AO  Al  A2  output_flag 

shoal.max.AT  26oct63  shoal  39.200n  118.380w 

piledr.max.AT  02jun66  piledriver  37.230n  116.060w 

rubis.max.AT  20oct63  rubis  24.000n  5.0(X)e 

saphir.max.AT  27feb65  saphir  24.060n  5.030e 

here  the  first  2  lines  give  label  (format  a80),  event-event  spatial  correlation  parameters  (free  format)  as 
well  as  the  output  terse  level  control.  The  inter-event  correlation  coefficient  is  defined  as  follows:  AO/2 
exp(  -(d/Al)'2  )  +  AO/2  exp(  •(d/A2)‘2  where  AO,  Al,  and  A2  are  input  parameters,  and  d  the  distance 
between  two  evnts  of  interest.  An  enupirically  determined  parameter  set  (A0,A1,A2)  » 
(0.8,20km,S00km)  has  been  used  in  most  of  late  calculations.  Each  of  the  remaining  lines  in  the  input 
file  specify  the  event  file  name,  date,  event  name,  and  geodetic  coordinate  of  the  event  with  format 
(al9,lx,a7,3x,al0,3x,f6.3,al,2x,f7.3,al). 

Each  event  file  contains  a  list  of  stations  as  well  as  the  corresponding  measurement  with  format 
x,f4,2)  as  showr;  in  the  followipj  sample  file: 

aae  >  5.36 
aku  5.53 
alq  5.62 
ant  <  5.39 
aqu  4  71 
are  <5.16 
atu  5.28 
bhp  <5.47 
lah  <  5.64 
Ipa  <  6.27 


Sun  Rclea.se  3.4  Last  change:  Sept  87  1 
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pda  <  5.93 
qui  >  6.23 

DIACNO.STICS 

Output  buffer  fort.38  contains  diagnostic  message  of  SVf)  routines  which  would  indicate  how  SVD  han¬ 
dles  matrix  factorization,  inversion,  and  synthesis. 

Also,  larger  output  flag  will  give  detailed  intermediate  output. 


SKK  Al„SO 

VSC-TR-82-12,  TGAL-TR-83-06,  and  TGAL-S6-01. 
•Source  codes  are  stored  in  "geotech/sre/emd/magy/GLM. 


ALTItORS 

R.  H.  Shumway,  R.  R.  Blandford,  and  R.  A.  Wagner  initialized  the  code  design  during  1982-1983. 
K.L.McLaughlin  revised  and  extensively  utilized  the  code  during  1983-1985.  R.-S.  Jih 

debugged/documented  the  code  and  implemented  SVD  routines. 


HUGS 

miglmd  needs  huge  memory  and  storage  to  execute  when  large  data  set  is  to  be  handled  and  inter¬ 
correlation  correction  is  desired.  While  mlglm7  has  been  used  lately  to  invert  300  events  and  150  sta¬ 
tions,  current  upper  limit  on  number  of  unknowns  (i.e.  stations  plus  events)  that  mlglm6  can  handle  is 
about  280. 
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